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ABSTRA CT 

Efforts related to the utilization of tunneling junctions for 
the conversion of radiative power to electrical power are reported. 

The theoretical foundations for these particular devices is presented 
along with a discussion of many of the practical problems associated with 
the implementation of such devices. 




I INTRODUCTION 

The enclosed report is suinbitted as an in depth review of work 
supported by N.A.S.A. grant number NAG 3-88. It in essence is the cul- 
mination of a number of years of investigation into the utilization of 
rectifying contacts for the purpose of converting either a monochromatic 
or broadband (solar) spectrum into electrical current with efficiencies 
approaching those achievable or exceeding low frequency classical 
rectifiers. Although such conversion efficiencies are not achievable in 
terms of presently available devices, the theoretical understanding and 
resulting predictions are beginning to indicate junction and coupling 
requirements which would have to be achieved to demonstrate the possibility. 
Recently also related areas of development, such as surface enhanced 
Raman generation, have stimulated interest in the coupling to micro- 
structures that should prove relevant, particularly with regards to the 
coupling of the incident radiation field to the junctiCT devices. The 
study of tunneling devices for conversion should also play a relevant 
role for p-n junction and other devices which have already demonstrated 
conversion efficiencies in excess of 10%. Such questions as the role 
of phonons and the influence of phonon collisions on the open circuit 
voltage problem might be illuminating from a fully quantum mechanical 
point of view as is necessary when tunneling is considered. Also since 
tunneling junctions have already been demonstrated to be useful with 
various types of cells, the question is relevant. 

The report is organized in a manner which presents first of all the 
theoretical foundations of some of the quantum mechanical aspects of the 
problem. These are particularly important for many of the concepts which 
wi'l be necessary to pursue. This section is the basis of an understanding 
of the relation between the quantum mechanical and classical picture of 
the rectifying process and in particular its relationship to the multi- 
photon theory of rectification. At optical frequencies this ties together 
the transition rate across the junctions and its relationship to the 
elastic current voltage characteristic. This is based upon a stationary 
state model of tunneling which provides the simplest way of tying tlie photon 
assisted tunneling theory of Tien and Gordon to the detailed interactions 
which occur in microjunction devices. Much of this effort is continuing 
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since the connections made are recent results. 

In the second section of the report we discuss the particular modes 
which will be of interest in operating devices. This involves the solution 
of dispersion relations for the inodes of the micro-structure junctions 
and in particular, the derivation of the loss coefficient associated 
with some of these modes. This is a particularly crucial problem since 
the modes turn out to be extremely lossy and in addition travel very 
slowly. An understanding of this dispersion curve for various frequency 
regions is important. 

Open surface (one metal and an overlaying dielectric) modes are 
also discussed and although they have more reasonable propagation lengths 
and loss coefficients, they still are highly lossy, 

In section IV of the report we discuss the coupling of open surface, 
free space and confined or gap modes. This involves basically three 
types of structures, the prism or grating, antenna and aperture, or, 
taper. Of these the first is very difficult because of the short 
periodicities required, The second is probably feasible at optical 
frequencies in spite of the short propagation lengths required. Taper 
coupling although difficult in terms of present day fabrication techniques, 
should become feasible and could eventually be incorporated into an antenna 
coupling scheme. 

The next two sections of the report pertain to experiments which 
have been undertaken to investigate some of the underlying principles 
and coupling techniques. This begins with Chapter TV, section 5 des- 
cribing a heterodyne diagnostic experiment. In Chapter V the fabrication 
of our junction devices which we are attempting to utilize for most of 
these experiments is described. 

We conclude the report with our perception of the prospects of obtaining 
practical devices which could display a relatively high conversion effi- 
ciency. We conclude that the theoretical possibility is present, however 
the fabrication techniques as well as the ability to obtain junctions 
exhibiting the desired characteristics are not yet at hand. Many of 
the ideas developed in the present report have been extended more recently 
and appear in the renewal proposal. The overall opinion is that whereas 
a few years ago it was uncertain as to how tunneling junctions were 
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managing to mix signals together, to couple the radiation field into the 
junction, and to exhibit the bias dependencies which they have, it is now 
felt that a rather thorough theoretical understanding of the physical 
principles Involved, the coupling, the propagation of the field along 
surfaces and into the junction, and the rectification characteristics 
of the junction devices has been obtained. The main problems requiring 
solution are those related to the fabrication of the structures required 
both for the tunneling and for the coupling in order to experimentally 
demonstrate that the theoretical modelling is correct. These problems 
are a challenge to microfabrication techniques which are presently 
evolving, such as molecular beam epitaxy, electron lithography, and 
chemical vapor deposition. 

The conclusion is that three aspects must emerge. First is the quantum 
mechanical basis of rectification theory and the inclusion of aspects 
such as phonon collisions and electron-electron collisions. Second, a 
junction displaying asymtiietry (no back current) must be theoretically 
analyzed and be experimentally fabricated. Lastly, coupling devices for 
the radiation field must be fabricated to demonstrate the possibility of 
significant deposition of energy into the region of the junction. 

While this report contains results pertaining to the direct conver- 
sion of the electromagnetic field to electron current, some of the basic 
findings have resulted from efforts associated with related investigations 
of the coupling between the free space radiation field and surface electro- 
magnetic waves and the coupling of the latter to electron tunneling 
currents. These sources have been NSF Grant No. ECS-7923877-02 pertaining 
to the investigation of tunneling devices for optical comriiunication, the 
Joint Services Electronics Program which has provided support for submill- 
meter and infrared mixing and detection and in particular the comparison 
of metal-barrier-metal tunneling devices with other types of junction 
structures. The Jet Propulsion Laboratories has provided support for 
surface wave amplification devices which utilizes the same basic expressions 
as for absorption. They have also provided some support for conversion 
work which has been described in a final report submitted and does not 
significantly overlap the presently reported effort. 
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With regard to notation, every attempt has been made to be as 
consistent and as clear as possible. The equation referencing has been 
made simple where possible. Where a single number is found reference is to 
an equation within the same section. If two numbers are given the first 
is to the section number of the same chapter and the second the equation 
number. All three numberfj are given when reference is made to a differen 
chapter. In section (I'l/.J), the section number was deleted. 

In final printing it was noticed that in some cases although every 
attempt was made to retain a consistent notation throughout, there is in 
some cases double notation. The most important of these we list below: 

1) 3 or k is used to represent the propagation vector component 
in the direction of propagation of the surface electromagnetic 
mode . 

2) 3=3' ‘ 13" = k = k' - ik" 

3) Y es well as v is used to represent the propagation vector 
component in the direction transverse to that of the propagation 
of the mode. 

4) S or 6 is used to represent the relaxation time of the free 
electron cloud. 

5) Barrier thickness is given by s, d, or b. 

6) The elastic junction current is given by <(V0)) in Sec. 11,2. 

I 1 in Sec. 1 1. 3 and simply I in Sec. 11.5. 

e I i k 

7) References are indicated by [] or 
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II THEORETICAL CONSIDERATIONS FOR RECTIFICATION 


II-l Introduction 

For the application of tunneling devices for the purposes of converting 
between modes of the radiation field and low frequency modes of lumped 
circuit devices, a suitable model describing the tunneling process is 
required* The model must include the effects of photon assisted tunneling 
as a classical analysis does not apply when the photon energy is large on 
the scale of the nonlinearity. The current theory, known as the transfer 
Hamiltonian model, has the dual disadvantage of being inapplicable in the 
thin barrier limit, at least for practical considerations, and an inability 
to treat barrier height modulation produced by the radiation field. A 
theoretical approach of Tien and Gordon, and used by Tucker relates the 
inelastic current resulting from an interaction of the radiation field 
to the low frequency elastic current. This relationship will become 
central for calculations of the potential efficiencies available. However, 
it was deduced rpl by considering the detailed interaction within the 
barrier region but by discussing the relative modulation of the wave 
functions on one side of the barrier with respect to those on the other 
side of the barrier. 

We have reviewed this problem and have been able to tie the results 
of Tien and Gordon to the interaction itself, This has given us an accurate 
picture Of the interaction process and in particular, its influence upon 
the barrier. We can thus include processes which do not necessarily re- 
sult in the Tien-Gordon formulas and which will become important in dis- 
cussing the conversion problem. 

This particular approach involves an extension of the stationary state 

model' " ■' to the case of an a.c. bias excitation, and does not suffer 
from the above limitations. 

1 1-2 The Stationary State Theory of Tunneling 

a) Static Bias Case 

In the stationary state approach, one initially finds the eigenstates 
of the Hamiltonian describing the complete structure. The exact form of 
the barrier potential is a complicated and generally unknown quantity, and 
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is a strong function of the material properties and interface quality. 

One generally has to assume a fonn for this potential, which represents 
a phenomenelogical input. It is true, however, that the tunneling pro- 
bability is relatively independent of the barrier shape, and one can 
generally assunie a phenomenological barrier height and width for an assumed 
rectangular barrier. A rectangular barrier allows straightforward deter- 
mination of the eigenstates, although closed form solutions exist for 
triangular, trapezoidal, parabolic and perhaps other simple geometries. 

For the treatment of field emission, for instance, a triangular barrier 
is most appropriate. Furthermore, for arbitrary barrier shapes with 
gently sloping interfaces, the WKB approximation can be employed. 

Initially this chapter will restrict discussion to tunneling between 
free electron metals. This is not essential and certainly other junction 
configurations will be important for conversion devices, The metals are 
characterized by the chemical potential p, which governs the density of 
electrons in such a manner as to preserve charge neutrality with the 
positive lattice in thermal equilibrium. When assembled into the tunneling 
structure, the Fermi levels must coincide in thermal equilibrium (this 
condition actually defines thermal equilibrium). The energy difference 
between the conduction band edges of the two metals is therefore well 
defined, as is the height of the effective barrier as measured, e.g. 
from one conduction band edge. The structure is shown schematically as 
an energy-displacement diagram in Fig. Il-la. An applied bias potential 
will shift the Fermi levels with respect to each other by an amount eVg. 

This assumption is due to the large effective screening in the electrodes, 
which implies that the electrons will "see" the field due to the applied 
potential only in the barrier region. This further assumes that there is 
appreciable reduction of the electron wavefunction amplitude in the barrier, 
i.e. the charge density and hence the screening are much lower there than 
in the electrodes. The energy diagram of a d.c. biased (conventionally 
positive) junction is Shown in Fig. Il-lb. From this model, one can 
inwediately write down the structure Hamiltonian from which the eigenstates 
are to be determined. 

If interactions between the electrons are neglected, the eigen- 
functions can be separated into two disjoint subsets, one associated with 
each electrode, and treated as independent thermal ensembles characterized 
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by the left and right Fermi energies. Therefore* once one assigns a 
suitable basis of structure eigenstates* the expectation of observables 
such as current density can be computedr 

In the ideally isolated metal electrodes, each energy level corresponds 
to a standing wave Bloch state occupied by an opposite spin pair of 
electrons. Ihe energy density of such states is 0 ^( 0 ) and Pp(c) in the 
left and right metal* where c is measured from the respective conduction 
band edge. 

For the real tunneling system, we wish to take into account the fact 
that for an electronic state which has a large amplitude in the left 
electrode, there is a finite albeit small amplitude for the state in the 
right electrode. It is assumed here that the eigenstates have the following 
formj the states which originally were standing wave states in the left 
electrode are now matched to a traveling wave in the right electrode carrying 
flux away from the barrier* and similarly the right electrode states are 
now matched to a traveling wave in the left electrode carrying flux away 
from the barrier. The states, of course, are functions of the applied 
bias voltage Vq across the barrier. We shall denote the new states corres- 
ponding to the original eigenvalues of the left electrode (1 c^(Vq)>}, 
and those corresponding to the eigenvalues of the right metal 
{|b^(Vo)>}. The {e^} and {cp}» measured from the respective conduction 
band edge, are distributed according to pj_(c) and Pp(c). 

The Hamiltonian as a function of coordinate x is 


H 



+ 4^(x) + eVQ 


f(K), 


(II. 2.1 


with 


<Hx) 

f(x) 


x<0 0<. x£b x>b 

6 4>jj(x) 0 

1 fb(x) 0 


ORiti'ff'ji**;, a i,, c -■ 

OF POOR QUALiTY 


(II. 2. ^ 


where 6 = and cp|^(x) is a complicated function representing the 
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effective barrier. The function fb(x) represents the normalized bias 
potential in the barrier. We will take the average total barrier height 
to be + cVq/2, i.e. the Hamiltonian from which the eigenstates will 
be defined as 

x<0 0<x<b X’>b 

2 eV 

^ ■ liii ^ + 4) 0 (J 1.2.3) 

Eigenenergy of H: e + oVq + 6 , for |C^(Vq)> and for !b^(VQ)) 

It should be reiterated that although a rectangular barrier is assumed 
other approximations, e.g, WKB, may be used. 

For high bias voltage Vq s tj), the barrier potential is more 
appropriately described as triangular or trapezoidal, hence the assump- 
tion that oVq < <p for these eigenstates. Structure energy levels are 
referenced (arbitrarily) to the conduction band edge of the right metal. 

The unnormalized (~) eigenstates have the form 


Eigenstate 

X < 0 

0 <, X b 

IW 

E 

C e"*^^+D 

E C 


T' e"^l^ 

' e+eVQ+6® ' 

E+eVQ+6 ^e+eVQ+6^ 


X > b 



ik.,(x-b) 
T e ** 

E 



-ik-(x-b) 

*3 . r% 1 

ik2(x-b) 


where 
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kj = /2meA»^ 

K = r^m((|.-eVQ/2-6-c)/h^ 

kg = /^m(c+eVQ+6)/h^ 


*s#* 


P6 POOR 


V 8 . 


iT 


and it may be assumed c € {e^} for Ic^(Vq)) states, and c + oVq + i> 


Matching amplitudes and derivatives at the boundaries x = 0 and 
X = b implies 


2 ik^(K-il^) 


K+ik- 

-2Kb/ 


-Kb 


R_ = 


(K+ik^)(K-ik3)-e"^''‘'(K-iki)(K+ik3) 


T, = 


4ikjKe 


-Kb 


where 


= -(K-iki)(K-ik3) + e'^'^‘^(K+ik^)(K+ik3). 

Transforming kj -»■ k 3 , *<3 k^, yields the corresponding quantities 


'^c+eVo+6- '^G+eVg+i have 

the important properties: 


^3(1- |f>E+eVg+«l^^ ' hj . 

and we observe also 

^ = ^e+eVQ+6 ’ 


f 

y 

(II. 2 . 5 ) 

Ur) 


(11.2.6) 

(11.2.7) 


(11. 2. 8) 

(11. 2. 9) 

(11. 2. 10) 
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so 

kj 

^c+eV(j+6 
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The states derived have appeared in the literature previously^^’^^ . 

A clear discussion of orthogonality, however, has not been oiven 
and is important for the following a.c. interaction considerations. One 
might assume that the wavefunctions are orthogonal outside of the barrier 
region, where they assume plane wave character, however this is inaccurate 
since the barrier region will always contribute an overlap. It might 
then be argued that the overlap would become insignificant as the normali- 
zation length tends toward infinity, as the overlap is independent of the 
normalization length and hence the contribution should decrease as 1/L. 
This argument does not provide an adequate degree of orthogonality, as 
the overlap is always of the same degree (1/L) as the current operator 
martix elements. It can, in fact, be shown that in the continuum limit 
the wavefunctions are exactly orthogonal. 

Consider for example orthogonality among the lc> states. Assuming 
an infinite structure and using 


e’^^dx = 7f6(A) + iPd), one has 

0 A 


< c^,(Vo)lc^(Vo)> = ,i(UlRj2)6(kj-ki) + n|Tj.l^ Mkj-k^) 


+ barrier contribution 


( 11 . 2 . 12 ) 


It is arduous but straightforward to verify that 

the principal part terms cancel the barrier contribution. 

3k, k. 

6(k3-k') = 6(ki-kJ)^ -F'^Ai-ki) 


Using 

one obtains 


< c^.(Vo)|c,(Vo)> = 27r6(kj-kj) 


(1 1. 2. 13) 


- 12 - 


ORIGINAl PA6u 

Similarly one may show OF POOR QUALrrV 

<b^.(Vo)lb^(Vo) = 27i6(k3-k') (II. 2. 14) 

and as T^tRg<S(kj-kJ) + R^*Tj,6(k3-kJ) = 0, one has 

<b^.(Vo)| c^(Vq)> = 0 (II. 2.15) 

Orthogonality holds exactly only in the continuum limit, but the 
discrete formalism developed in this chapter is applicable with the 
understanding that evaluation of expressions should be done in the L-^*' 
limit. Adopting the normalization <c^(Vq)| ■;;;^^c^(Vq) j , 

<b^t(Vo)| = which is usual, one may return to the summation 

notation with normalization 


<bj,,(Vo)|Cj.(Vo)> = 

and we have the associations 6(k^-kp = 6|^ |^, 


(11. 2. 16) 

(11. 2. 17) 


^6(k3-k3) = 6|^ 1^, as L-ho. In a practical sense, going to the continuum 

limit involves only changing summations to integrations in the usual 

manner, i.e. I p(£)dc, and taking the domain of the wavefunctions 

£ 

to be (-“,“). The change has no effect on the current operator matrix 
elements. 

Charge neutrality in equilibrium is maintained, I.e., the electron 


density in each metal is the intrinsic density at zero bias. For 


instance, in the isolated left metal for a given energy the electron 


density is n(k) = -^ . 
we have n(k) = i((l 
right metal as well . 


Forj^the left metal in the structure with Vq 
IT' ^ * This is true in the 

^3 ' C+(V ' II 71 



f 
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The current operator is 
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(II. 2. 18) 


where x' is now taken to be the coordinate variable upon which the wave 
functions depend and x is the position of observation. 


Evaluating the matrix elements, one finds 
<Cj,(Vo)|jlc^(Vo)> =|^(k3 + k3) • 

^ ^c+cVq+S ^ ^ I ^c+e Vq+iS ^ 

2m 1 2^^|P^|2 


(II. 2. 19) 


(II. 2. 20) 


The fact that the final momentum must be real precludes current flow 
below the conduction band edge of the electrode into which the electron 
tunnels. Multiplying the matrix elements by the respective standing 
wave state densities (including spin degeneracy, and where 0(c) is the 
Heaviside function) 


pl(e) =%Tf-e(.) . 

HTT 1 

pp(c+eVQ+6) e(c+eVQ+6) 


one has immediately 


PL(E)<Ce(Vo)|j|c^(Vo)> 


(II. 2. 21) 


(II. 2. 22) 


=0, (II.2.23) 

the condition of detail balance. This means that for any structure 
energy level, if both 1 Cj.(Vq)) and equal occupation 

probabilities., then there is no net current across the barrier. This 
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ensures that there is no net current flow at zero bias. 

In the second quantized form, ORIGINAL *‘3 
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H = I ^ (II. 2. 24) 

•L 

where the creation operators c^, create |c^(Vq)> and lb^(VQ)> eigen- 
states with eigenvalues e + eVg + 6 and e, respectively, i.e., 

HlC|.(Vp)> = (e+eVp+«)ri:^(Vo» (11.2.25) 

H|b^(VQ)) =Elb^(VQ». (II. 2. 26) 

One generally utilizes the complete set of eigenstates to form unperturbed 
Green's functions, as this facilitates the treatment of, e.g., phonon 
interactions [4]. Here all interactions are neglected, thus there is no 
need for the introduction of Green's functions at this point. 

The time dependence of the field operator is written explicitly 

c Ic (V„)> e 
'^L 

+ I bJb^(V(,))e'''='^'” (II. 2. 27) 

so that ij; satisfies ihi|j = (Heisenberg picture) and the creation 

operators are time invariant. 

It should be pointed out that here Dirac notation is used 
to represent the Schrodinger wave function. 

Closure of the Dirac brackets implies an integration over the coordinate 
variable (x') upon which the wave functions depend. The density operator 
is 

P = I e“^^’ (II. 2. 28) 


-15 


where 
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r'n.j-, 


!;-r 


K = H- i (MR + eVo)<c, - I PR^'b,. 
and 

Z = Tr e"®'^ . 3 = A 

Kgl 

Thus 

K = J (^-\)4,ve * I • 

"L =^R 


( 11 . 2 , 29 ) 


(11.2,30) 


( 1 1 . 2 . 31 ) 


This expression for K shows that the bias dependence is in the states, 
i.e., the density operator is form invariant with changing bias. The 
second quantized current operator is J = or 


J= Z , c+c^,<c^(Vo)U|c^.(Vo))e-’(^'-^)‘/'' 
R’ R 


+ 


+ 


I ,^A'+eVn+6^^c^'^0^ 




_J ^c+eVp+6 ^c' ^ ® 


-i (c'-c)Vh 

(II. 2.32) 


Spin indices are suppressed for notational simplicity, as the 
Hamiltonian and density matrix are trivially spin diagonal. The current 
operator matrix elements are independent of spin, hence the summations 
over c,g' imply additional summations over independent spin variables. 

In this chapter all spin summations can be performed by including spin 
degeneracy in the density of states, thus one may adopt degenerate 
state densities and neglect spin throughout. 

For the static bias case at hand, the current operator may be evaluated 
anywhere in the structure, and the total current J(Vq) = Tr(pJ) may be 
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Ininedlately written down. 


<J(Vq)> = I 




oRiGirifU rA :% : 
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' <c,(Vo)Mlc^(Vo)> 


X 1+e 


-3(c-p|^) 


He 


(II. 2. 33) 


or in the continuum limit 


<J(Vq)> = dcpj_(c)f(r.-(j|_)<Cj,(VQ)lj|c_(Vg)> 


0OO 

dpp^(e)f(c-p„) <b^(Vj|Ho|b^(VQ)> (11.2.34) 

J mtOO 


where f(x) =- ^ is the Fermi function. Detail balance (Eq. II. 2. 23) 

1+g-px 

may be used to combine the two integrals. Then, 


mOO 

<J(Vq)> = dGp^_(e)(f(c-y|_)-f(c+e\/Q-M|_)) < c^(Vq) Ijlc^(Vg)) , 

(II. 2. 35) 

recalling that - Ml " 

Given p. (c) = — ^ 0(e), 

^ Trkjh^ 


and evaluating 


<c,(Vo)|JlcyVo)> 


4aikjK^(l<3 + k3)e'^'^^ 
m LlF 


one obtains 


8e r kjK^(k 3 + k 3 )e“^^^ 

<J{Vg)> - 3 g 3 (f(,-yJ.f(c+eVg-yL)). 

'''c' (II. 2. 36) 
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This agrees with the result derived by Harrison [10] to lowest order g 

in the tunneling exponential Fluctuations in the tunneling current^^^^ f 

can be treated similarly^^^ and give the results of Ref, (12) to a good 
approximation. 

b) Nonstatic Bias Excitation 

Previous to Ref. (7), the stationary state approach was restricted to 
static problems, i.e., where the applied bias was time independent. 

Although scattering from time dependent impurities was considered in 
a stationary state approach, Ref, (5), the response of tunnel junctions in 
the presence of a radiation field (effective nonstatic bias) was accom- 
plished only in the transfer Hamiltonian format^^^*^'^’^^^ . In this 
section, the stationary state method of Ref. (7) for photon assisted tunneling 
will be described. Due to the quantum nature of the radiation field, 
expressions for the tunneling current revert to summations over the 
number of quanta absorbed or emitted into the field, and nonclassical 
behavior is exhibited. By "classical behavior," we mean the predicted 
response of the junction when Vq V(t) in our previously derived ex- 
pressions in the static bias case. Throughout the present section, 
comparison will be made to results obtained from the transfer Hamiltonian 
format, and certain discrepancies will be noted. The discrepancies are 
the result of the change in barrier height due to the applied field being 
neglected in previous analyses, as will be described. The ability of 
the stationary state model to predict the effects of photon assisted 
tunneling demonstrates the complete utility of this model for independent 
tunneling charge transport. 

The approach described here is to find time varying eigenstates for 

I 

a structure with a sinusoidal bias component. The current of the time 
varying eigenstates is then summed in the usual manner, yielding the 
current response. The source of the time varying bias is not specified, 
and therefore could be due to a radiation field or plasmon field, within 
a certain frequency limit to be described, it is assumed that the fre- 
quency is low enough so that the density operator is time invariant, which 
is equivalent to assuming the maintenance of thermal equilibrium in the 
population of the {|c>} and (jb >} states. We are neglecting bias. field 
penetration into the electrodes so that c, the kinetic energy of the 
incident electron, is still a valid label. Approaching visible frequencies, 
the assumption of rigid motion of the Fermi levels (rigid occupancy 
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hypothesis (4)) breaks down, as does the assumption of a piecewise 
linear potential taken for V(x,t) « V(t)f(x). One should use perturbation 
methods for excitation at these frequencies, but the expressions observed 
here are still valid to lowest order in the a.c. field amplitude. The 
density operator therefore retains its form.and time dependence is intro- 
duced through explicitly time varying states. Consider the |Cj,(Vq)) 
states, which satisfy H|Cg(VQ)) = (c+eVg+a) |c^(Vp)) for static bias Vq. 

For an explicitly time varying bias voltage V = Vq + cos wt, |c^(V)> 
will no longer be associated with a simple exponential time dependence, 
yet if one now includes the time dependence in the state one must have 

H|c^(V(t))> = ih |c^(V(t))) . (II. 2. 37) 


Thus, states lc^(V(t))> must be found such that 



+ 4>(x)+e(VQ+V^ cos djt)f(x))|c^(V(t))> 


= |c^(V(t)» , 

(II. 2. 38) 


If ]c (V(t)) and the corresponding |b (V(t)> can be found, the a.c. 
current response neglecting field penetration may be computed without ap- 
proximation. To date, exact solutions have not been determined, but an ap- 
proximate solution may be found as follows. One assumes a solution to Eq. 
(2.38) of the form 


|c^(V(t))> 



-i (t:+eVg+iS+nt.ii)t/I' 


9 


(II. 2. 39) 


where the |c[!>are normalized, time invariant wave functions not necessarily 
orthogonal on n. Substitution into Eq. (2.38) above implies, because of this time 
invariance, that 


( ^ + 4>(x) + eVr,f(x) - (c+eVn+5+nrio)jA„ |c[l > 


n' E 


eV. 




n-1 


(II. 2. 40) 
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2 

( ^ * (j>(x)’*-(eVQ+nli[d)f(x))ICp> * (i+eVQ+6+nliw) |c” ) 


then one can see that jc^ ) is simply |c^(VQ+nliw/e)> , corresponding to 
the static bias case but with bias VQ+nWe, The'i one must have 


•f(x)nliwA^|Cg(VQ+nhw/e)> 


eV. 


+ ’-p f(x)(A^^jlc^(VQ+(n+l)iiu»/e> +A^_^ lCp{VQ+(n-l}lisu/e)) ) - 0. 

(II. 2. 41) 


Where f(x)=l (the left electrode, f(x) vanishes in the right electrode) 


j ik X ik .x 

Ic (V«+nhw/e)> has the form — (e +R" e ), where the reflection 

' c u ^ t 

coefficient is a slowly varying function of n when 

Neglecting this dependence in (2.41), i.e. assuming the above inequality 

implies, 

eV, 


•"’'t, + i hV <Vi+Vi> = ° 


eV 

^n ^ ^n ’ 


(11.2.42) 


J_ being the ordinary Bessel function. The time varying ]c (V(t))> 
n ^ 

states are thus taken to be 

eV^ -i (e+eVA+nhui+i^)t/li 

(II. 2. 43) 


i.e., a weighted sum over the previously derived states for the static 
bias case, where the bias potential differs by some integer number of 
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pho'cdn energies. Recall that |Cj^(VQ+rt»w/e)) satisfies 


( 1“ +(j)(x)+(e\/Q4n|iuj)f(x))|c^{VQ+nl»w/e)) * (c+eVQ+nliw+6) jc^ (VQ+nl»w/e)> 

(11,2.44) 


Similarly, one can find time varying jb(V(t))> states. The requirement 
that the transmitted part of the eigenstates be projected out of the 
recursion relations, so that one may use the approximate invariance of 
the reflection coefficient on n as above, implies (with hindsight) that 
a slightly different forin must be adopted for lb, (V(t))> , A suitable 
form is 


b^(V(t))> 



■t eV 

™ cos wt'dt' 

t»4) 


I BJb") 


(11.2,45) 


I 


and upon substituting this into the appropriate Hamiltonian 
expression one obtains the recursion relation 

2 

( +(j>(x)+eVQf(x)-c-nfiui)B^!bJ > 
eV 

+ (II, 2. 46) 

If one requires 

( lijf +'l'(x)t(eVQ-n/u .•f(x))lb^> ^clbj;) , 

then 

eV 

nMf(*)-l)B„|b") + -/(f(x)-l)(B„^j|b"'^b+B„.j|b"'^> ) = 0. 

(II. 2. 47) 

Where f(x)-l = -1 (the right electrode, f(x)"l vanishes in the left elec- 
trode) the states are approximately invariant on n, as it is evident that 
|b^ >>= |b^(VQ-nhw/e) > , thus implying that in this limit one should take 




eV. 




» J „ Thus, we adopt 

n -n ' hu> 


b_(V(t))> * e 


■1 


J -jpcoso.fdt' p^(^)lbJVo+nWe)>c'*<‘-"'‘"’>‘/'' • 


(II. 2. 48) 


Note that <Cg(V(t)) lCg(V(t))) =(b^(V(t}) lb^(V(t))) « 1, so that as stated 

previously the density operator retains its tinie independent form, however 
now the creation operators create the explicitly time varying states. 

It should be stressed that our solutions are approximate, and are valid 
only in the limit hw<< t|i, pj_, This inequality will be assumed 
for the remainder of this section. 

Our concern here is with the current, so the relevant (diagonal) 
part of the current operator is 

. 1 . eV« eV„ 

j = I cic I J (-f-C)J 

J: c £_^_, nMiuj' n'Miw' 

ilfll 

<c^(VQ+rihw/e)ijlc^(VQ+nhw/e) ) 

J. eV„ eV^ 

i ^ £r.^n t n' hoi ' n ' hw ' 

C|^ n,n' 

<bJVn+nhw/e)|jlb (V^+nhw/e) > (11.2.49) 


Evaluation of the trace y 'Ids in the same manner as for Eq. (2.34) letting 


n' = n+m 


eV eV 

e^,n,m 


< c^(Vo->-(n+m)ha)/e)|jjc^ (V^+nhw/e) ) 






E|^ , n ,m 


< b^(VQ+riiuj/e) 1 j |b^(Vo+(n+m)tuo/e) > e 


-in)o)t 


(II. 2. 50) 


The a.c. current is not a conserved quantity, (i.e. <()> f 0), 
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implying the existence of oscillating charge density. For instance, 
the |c^> matrix elements have sinusoidal spatial modulation in the 
right electrode. This is due to modulation of the electron wavefunctions 
on passage through the time dependent potential. For small w where 
classical circuit concepts apply, this spatial variation becomes arbi- 
trarily small. 

For the m s 0 (d.c.) terms, one can again use detail balance, i.e., 


P|_(c) <c^(VQtnfia)/e) 1 j |c^(VQ+nl»u)/e)> 

+p^(e+eVg+nlwt6) < ^+eVQ+nlw+S<"o+"'‘“''«) I j l‘>c+eVc,.nhu,+f < 


= 0 


(II. 2. 51) 


to combine the two integrals. One obtains 


-00 


' 5! I p(_(e)dc(f(c-pL)-f(c+eVQ+nha)-ML)) 


eV 

'^n ^ c^(VQ+nbo)/e) I j |c^,(VQ+nlWe) ) , (II. 2. 52) 


or 


0 eV„ 

I 4 <J(Vo+iir,a,/e) ) 
n 


(II. 2. 53; 


where <il(VQ+nhw/e)) is the current for static bias potential VQ+nbw/e. 

This is identical to Tucker's expression for the d.c. quasiparticle 

tunneling current derived from the first order transfer Hamiltonian 
( I5j 

model' ■ . It predicts the famous steps in the I-V characteristic due 

(13) 

to photon assisted tunneling ' . 

Eq. (1 1. 2. 53) shows that the index n can be interpreted as the number 
of photons absorbed in the structure and that the current is a superposition 
over the possible multi-photon emission and absorption processes. The 
total absorption rate y is thus given by the expression 

9 eV„ 

Y - I J(VQ+nlio)/e) /e 


(II. 2. 54) 


One can obtain expressions for the a.c. terms as well. This has 
been done as well as a consideration of the Kramers-Kronig relationship 
between the imaginery and real parts. Since these are not immediately 
relevant to the conversion problem they will not be discussed here. 

The primary results of the present section are the two above 
formulae II. 2. 53 and 1 1. 2. 54. They tie together the absorption aspects 
of the problem with the current voltage characteristics. In recent work 
we have been utilizing these expressions to investigate the conversion 
efficiency available from a p-n junction as a function of frequency. 

The tie in to the transfer Hamiltonian approach is extremely impor- 
tant since It allows one to begin with a microscopic interaction potential 
and discuss absorption and emission phenomena within the junction region, 
and within the limits of validity to obtain the Tien-Gordon results, which 
are obtained by consideration of the interaction outside of the barrier 
through a phase modulating gauge term. 

There are two aspects which allow this to occur. Firstly Eqs. 

(2.43) and (2.48) show that the n photon absorption and emission processus 
occur with an effective barrier height equivalent to a shifting of the 
d.c. potential to Vq +(-) nhw/e for the c states and secondly Vq -( + )nli(i)/e 
for the b states. 

This is what one would expect from an intuitive view as shown in Fig. 
(II-2) for the bias convention. Considering Eq. (11,2.5), from a transition 
point of view one can picture the final state as a b-state which for 
emission has an energy lower than the initial c-state by nliw and one 
might expect the K to differ for the two states. However for the former 
c-state the value of Vq must be decreased by nhw/e. Since the electron 
is undergoing emission from that state the effective barrier increased 
by nhw/2. On the other hand the final b state is first of all lower 
in energy by nhw and absorption occurs from it, thus for K the effective 
barrier must be increased by nhw/e. If we use these two facts for the 
final b-state in Eq. (2.5), we observe that the K value for the b-state 
is equal to that for the c-state, and is equivalent to a net increase 
in eVg by nh£jj/2e as observed above. 

These aspects make the matrix element which one could use in the 
transfer Hamiltonian approach different from that conventionally used 
for photon interactions. However it will be seen that the integration 
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over the densities of states to deduce current-voltage charateristics of 
various devices for practical consideration will become simpler. It is also 
this aspect which allows one to relate absorption to the current-voltage 
characteristic. 


II-3 Model of Absorption and Emission in Junctions to First Or* 
in Perturbation 

We now establish directly that to first order in perturbation the 
net rate of stimulated transitions into confined structure modes of tun- 
neling junctions is proportional to a finite difference on the elastic 
d.c. tunneling current- voltage characteristic. In the low frequency 
(classical) limit, this reduces to the differential conductance. Since we 
have specifically done this for a metal-semiconductor junction, we shall 
considerthis case. The result and approach are however general. The 
approach demonstrates the relationship of absorption to the current- 
voltage characteristic (as discussed above) through the interaction 
within the barrier and thus provides insight into the relationship 

between the stationary state and transfer Hamiltonian approaches and 
of these to photon assisted tunneling* 


Considering Fig. II-2, electrons in tunneling from left to right can 
emit a plasmon during the tunneling process and terminate in a lower energy 
state in the right electrode. Alternatively, the electron can absorb a 
plasmon from the field and enter an excited state in the right electrode. 

To deduce the relationship between the absorption and emission 
processes and the elastic current we represent the plasmon interaction by 
a time varying potential eV^ cos (wt) across the barrier region, and 
neglect the interaction within the electrodes. A one dimensional struc- 
ture is analyzed, and the plasmon propagation constant is assumed to be 
negligible. A rectangular barrier is taken, and the structure v/ave 
functions are approximated to lowest order in the tunneling exponential. 

One can express the net plasmon emission rate to first order of 
perturbation theory as 



J .oo 


dE PL(E){p^(E-ho)|M^^^_,.^l^f{E-y) - f {E+eV^-hm-y) ) 


-P|^(E+ha)) lM^^^^l^^l^(f(E-y) - f (E+eV^ +hw - y ) ) } 

(II. 3.1) 
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where m is the chemical potential of the left electrode* and for para- 
bolic bands the state densities are given by 


I 

f 


Pl(E) = Re-a^ , Pp(E') = Re-S^- 
h h 


(II. 3. 2) 


If the energy E is referenced to the co'.iduction band edge of the left 
electrode, and we specialize to an M-I-S structure where the final state 
lies in the valence band, one has 

* ’ ''•f(E+eVo) = (ii+*-eV(,-E')y , (11.3,3) 

where a single mass parameter has been adopted for simplicity. The 
parameter A represents the energy difference between the semiconductor 
Fermi level and valence band edge and is positive for the degenerate 
semiconductor specifically assumed for Fig. 2. 

The stationary state model of the last section shows that the decay 
constants of the initial and final electronic wave functions in the barrier 
are equal and given by 

K = (♦-(E+E'+eVg)/2)''* 

This relationship will be seen to be important in establishing the 
desired connection to the elastic current expression. 

The matrix elements where 

■ - e 


(II. 3. 4) 


«E'E 


dx 


'^L,E' 1" m^ ■ ^"2m *^R,E , 


(II. 3. 5) 


can be expressed for a guage choice of A = E/iw as 

eV. 


Mr. 


r ^-Kd 


E*,E 2w ’L'R m 


(II. 3. 6) 


where 


-2ik. 


-2ik^ 

^R ^ imr: 


(II. 3. 7) 
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The barrier height <p is measured from the left electrode (metal) conduc- 
tion band edge. 

Equation (3.1) is directly related to the d.c. current in the absence 
of inelastic tunneling, which to lowest order in the tunneling exponential 
is 


where 



(V 1 ^ 


dEpL(E)pR(E)lM^^£|^f(E-p) - f(E+eVQ-p) 

(II. 3. 8) 


M 


E'.E 


- 

" 2iii 


►00 

i -00 



(VS(x-|) 


( 11 . 3 . 9 ) 


Upon examination of equations (1) to (7), it is apparent that the energy 
shift of the final state is equivalent to a change in bias potential, 
-hw/e for the first (emissive) term and +1)03/6 for the second (absorptive) 
term. Comparison of (1) with equation (8) for the d.c. elastic current 
yields the relationship 

which in the classical (o3 0) limit is 


/htov /Vr\2 


9 

9V 


0 


>el (''o> 


(II. 3. 10) 


This is precisely the relation one obtains from the theory of photon 
assisted tunneling ^^’^^^in the limit eV^«hu3 . The present derivation 
depends specifically upon equation (3.4) from the stationary state theory 
and is obtained by considering the interaction within the barrier. It 
is a specific example of the general validity of the Tien-Gordon equation 
when barrier height modulation is included in the theory as in the 
stationary state approach [ 7 ], as opposed to the simple phase modulation 
imposed on the left electrode wave functions in the transfer Hamiltonian 
model^^^^ . 

Equation (3.4) has a simple physical interpretation which can be 
obtained by considering in particular a transition from a left state to 
a right state hw greater in energy. Viewed as a stationary process, the 
electron receives additional kinetic energy in making the transition 
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to the right electrode. Thus the barrier height it experiences emerging 
from the right side of the structure is reduced by Tiw and the average 
barrier height is reduced by fiw/2. Therefore to be correct, modified 
wavenumbers in the barrier region corresponding to an effective barrier 
height 4)-fU)/2 for the left states participating in absorption processes 
and an effective barrier height (it+hto/Z for the corresponding right states 
should be used. For emissive processes the signs are reversed. These 
corrections are generally small but are important in establishing the 
equality implied in equation (9) and its relation to photon assisted 
tunneling. 

II-4 Relationship of the Stationary State Model to the Transfer 
Hamiltonian Model 

In the present section of this chapter, it will be shown that 
the expression for d.c. elastic tunneling current itself can be written 
in the form of a transition rate, using a relation between the structure 
eigenstates which will be derived (17). This correspondence relates the 
stationary state model to models based on perturbation theory, in parti- 
cular the transfer Hamiltonian model. 

It is assumed that the tunneling particles are independent and have 
free particle dispersion characteristics throughout the structure. The 
extension to the effective mass approximation is straightforward, and 
does not change the nature of the result, so for simplicity a single mass 
parameter is assumed here. Furthermore a one dimensional structure is 
assumed, as the extension to three dimensions is also straightforward and 
again does not alter fundamentally the result of this section. Given that 
a static bias potential Vq appears across the structure, an orthogonal 
set of time independent structure eigenstates are assumed to be known. 

As in section one, the 1 c^(Vq)> are eigenstates associated with the 
left electrode and carry flux to the right, and deep in the right elec- 
trode the lCg.(VQ)> assume a plane wave form. Similarly the |b^{VQ)> are 
eigenstates associated with the right electrode and carry flux to the 
left. The parameter c is the particle incident kinetic energy. It is 
assumed that the barrier potential is of finite extent, and outside of 
this region the states have plane wave behavior with incident, reflected 
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and transmitted components. The electrodes have equal thickness which 
will be taken to approach infinity, hence the actual behavior of the states 
in the barrier region is unimportant with respect to orthogonality and 
normalization. The normalization used here will be that used previously, 


< c^i(Vq) 1c^(Vq)> 


= 6 


ki'ki 


9 


where one may associate, e.g., r<S(k,-k{) -(-► 6,. as we assume L > 
Outside of the barrier region the wavefunctionS have the form 


X » 0 

ik„x 
Te ^ 

-ik«x ik^x 
e +R'e (II. 4.1) 

where again the incident kinetic energy subscripts have been suppressed. 

The states are associated with the standing wave state energy density of 
the respective electrode, which are, neglecting spin degeneracy, 


X <•< 0 

ik,x -ik,x 
/?T|c(Vq)) e ^ +Re ^ 

/2Tlb(VQ)) re 


|c^(Vo)> ^ Pl = 

it 1 i 1 




n - JL, 

'' “ TTh^ ^ 


(11. 4. 2) 

(11. 4. 3) 


For any two degenerate eigenstates of the structure Hamiltonian, it is 
easily verified that the current operator matrix element between the two 
wavefunctions is constant throughout the configuration space of the Ham- 
iltonian. Therefore one has <c(Vq) [ j |c(Vq)> , <b(VQ) ] j lb(VQ)> and 

<c(Vq) I j jb(VQ)) constant throughout the structure for states corresponding 
to the same structure energy measured, e.g., from the right electrode con- 
duction band edge (the energy subscripts have been suppressed). One may 
also consider the states jc (Vq)> , |b (Vq)> which are the complex conju- 
gates of the (time independent) 1 c(Vq)> ,jb(VQ)) . The complex conjugates 
are orthogonal to each other, but not to |c(Vq)), Ib(VQ)) with which. they 
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are linearly dependent. With these states one can prove the following 
relationships between the exact structure eigenstates; 


< c(Vp)|jlc(V(,) > 

= - < c*(V())|j|c*(VQ)> , 

(II. 4. 4) 

< c(V(,)|j|c(llQ) ) 

|( c (Vg)j,)lb(VQ)) I^Pp . 

(II. 4. 5) 

(b(V£,)|olb(Vo)> 

= -<b*(Vg)|j|b*(Vg)> , 

(II. 4. 6) 

<b(Vo)|j|b(Vo)> 

= - ^kc*(Vo)|j,;b(Vg)) |\ , 

(II. 4. 7) 


where j is the current operator and again all states correspond to the 
same structure energy. The proof may be obtained through straightforward 
evaluation of the matrix elements. As all of the matrix elements are in- 
dependent of position, one may evaluate < c(Vq) ] j jc(VQ)) , <c (VQ)|j|c (Vq)) 
and < c (Vq) j j jb{VQ)> deep in the right electrode, yielding 

1 ehk^ « 

(c(Vo)ljlc(V(,)> (II. 4.8) 

it it 1 eliko n 

<c (Vq)|j|c (Vq) > (II. 4. 9) 

<c (Vq) j b(VQ) > =--^(-ir)T. (II. 4. 10) 

(4,5) then follow immediately. Similarly, evaluating < b(VQ) | j Ib(VQ) > , 

< b*(VQ) j j |b*(VQ)> and < c*(Vq) I j |b(VQ)> deep in the left electrode leads 
directly to (6,7) above. The expressions that are most interesting for the 
purposes of this section are (5,7), as they relate the current of an 
eigenstate to an expression in the form of a transition rate, with j/e 
being the "interaction Hamiltonian." 

The connection between the first order transfer Hamiltonian and 
stationary state models is established mathematically through (5,7). 

In the transfer Hamiltonian formalism* one generally adopts tunneling 
matrix elements as derived by Bardeen [6] so that the model reproduces 
the result of the stationary state approach to lowest order in the tun- 
neling exponential. The prescription applied is to first find eigenstates 
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of the noninteracting L and R subsystems which are found by taking the 
barrier to have semi-infinite width, i.e. the L subspace consists of the 
left electrode and barrier for x ^ 0, and the right electrode and barrier 
for X < b constitute the R space* The eigenstates, therefore, are standing 
waves. The expression for d.c. tunneling current is an occupation weighted 
sum of the transition rates between L and R states, where bj/e, evaluated 
in the barrier region, is taken as the interaction Hamiltonian. 

Although this procedure does reproduce the d.c., current expression 
of the stationary state model to lowest order in the tunneling exponen- 
tial, the accuracy of higher order terms in the perturbation expansion 
is questionable due to the first order eigenstates employed. The eigen- 
states of the L and R subspaces of the conventional transfer Hamiltonian 
approach are approximations to lowest order in the tunneling exponential 
of the |c^(Vq)> , |b(VQ)) states (or equivalently, |c(Vq)> and |b*(VQ) )) 

above. Thus the fact that this procedure approximates the stationary 
state result is a consequence of Eq. (5,7). An exact correspondence is ob- 

lir lAr 

tained if one adopts [c (Vq)> states for L and |b (Vq)> states for R in 
the first order transfer Hamiltonian model. This gives: 

de P|_(c)(:^^(e+eVQ+6)(f(e-i.iL) - 





(II. 4. 11) 


which using (5) is equal to 
<J> = 


dc PL(c)(f(c-PL) - f(cfe\lQ-]x^_)) < c^(Vq) j j Ic^(Vq)) 

(II. 4. 12) 


where 6 is the difference in electrode chemical potentials This 

is the exact expression for the current from the stationary state 
model [?]. Replacement of Bardeen's approximate L and R wavefunctions 
with the exact ]c (Vq)) , Ib(Vg)) thus makes the two models equivalent 
with respect to the computation of the d.c. current. 

One may view the first order independent particle Hamiltonian 
model as a formal approximation of the stationary state model where the 


I 

f 
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use of hj/e as the interaction Hamiltonian is justified by Eq. (5,7), 

With the association of the first order and full structure eigenstates 
and Utilizing (5,7) one can regain the stationary state d.c, current 
expression from that of the first order transfer Hamiltonian model. This 
fact may be useful in extending the first order model to, e.g., extremely 
thin barriers. 

In this section we have described how the stationary state and first 
order transfer Hamiltonian tunneling models are related, in the case of 
static bias and independent particles. The relationship is a consequence 
of the properties of the structure eigenstates, in particular that the 
flux of an eigenstate may be expressed in the form of a transition rate. 
Using the relationship, one may construct an exact transfer Hamiltonian 
model for static bias and neglecting collisions simply by using the appro- 
priate structure wavefunctions instead of the commonly used approximate L 
and R wavefunctions. 

II-5 Theoretical Considerations for Practical Converters 

The basic relutionships developed in the previous sections which are 
most important as far as mode conversion in tunneling junctions is con- 
cerned are first of all that providing the d.c. current, inelastic (induced) 
and elastic (back-flow). Secondly one has the expression providing the 
transition rate due to the presence of the mode within the junction itself. 
These two are: 

(11.5.1) 

(11. 5. 2) 

One observes that to maximize energy absorption from the field the current- 
voltage relationship should satisfy I(VQ-fnriw/e) » I(VQ-nliw/e) for n 0. 
This is similar to the condition one would expect for classical rectifi- 
cation, although due to the quantum nature of the radiation the change 
in slope need not be abrupt. Only points on the current-voltage charac- 
teristic corresponding to integral numbers of photon potentials contribute, 
which at optical frequencies are separated by the order of one volt. For 


? eV 

^ " ^'^n I(VQ + nho)/e)/e , 
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the current-voltage relationship exhibited by an Al-AlgOg-Ai field 
emission junction, the rectification efficiency increases with increasing 
frequency. At visible frequencies one would expect high rectification 
efficiency as the current change can be several orders of magnitude over 
a photon potential. 

In order for energy to be extracted from the field, the hot electrons 
produced by photon absorption must be prevented from decaying thermal Ty 
within the electrode. Electrode band structure can be used to forbid 
decay through thermal izati on, for example if the fin^l tunneling state 
is a narrow impurity band in a semiconductor band gap, the energy loss will 
be less than the impurity band width. Such a configuration is shown in 
Fig. 11-3 (photon "pump") and is to be considered further. 

In order to extract net energy from the field it is necessary to 
compare the absorption rate given by the above formulae with the losses 
inherent in the system. The sources of loss can be divided into two types; 
those caused by the inefficiency inherent within the process itself, and 
secondly losses associated with other processes accompanying the direct 
conversion process. An example of the former is the harmonic conversion 
in the rectifier converter, since it arises as a result of the conversion 
process itself. The efficiency of such processes are given by the ratio 
of the useful transition rate over the total transition rate. For the; 
rectifier this is the ratio of the transition rate providing the d.c. 
component to the total transition rate. This aspect is dealt with in 
more detail in the recently submitted proposal. 

An extremely Important example of the second type of loss and one 
which must be carefully considered for the tunneling converter is the 
propagational loss associated with the particular electromagnetic mode 
which is inducing the transitions. This particular loss is given by 
the imaginary part of the propagation constant associated with the pro- 
pagation of the mode. This loss is discussed more fully in Chapter III 
pertaining to the structure modes of Interest. 

In order to compare the propagational loss with the useful transition 
rate one must introduce the loss coefficient for the conversion process. 
This is given by 


a = hwR/P 


( 11 . 5 . 3 ) 



where R is the absorption rate per unit length of the device and P is 
the power propagating in the mode inducing the transitions. The former 
is given by y/^ 2 * ^ interaction length of the converter. 

P the power is given in terms of the Poynting vector 


P « i. 


E X R)2 ds . 


To lowest order in perturbation y is given by Eq. (II. 3,1) which 
for the three dimensional case can be written as 


Y * 


« 2ir 


- f(E,+eV„-p) X 


(6(E,-Ef+hu) - S(E,-£f-l.u))) 


eV. 




h W\ 


fi uu 


e'2fjw 


(II. 5. 4) 
(II. 5.5:) 


the latter by Eq. (1 1.3. 9) 

One observes from Eq. (5.3) that for appreciable absorption to occur, 
the ratio of R/P must be as large as possible. Since R in the region of 
interest is proportional to the power propagating within the barrier, 
this implies that the greater the ratio of power propagating within the 
tunneling barrier to the total power propagating within the mode, the 
greater the value of a, and hence the greater the overall conversion 
efficiency possible. In a region of normal dispersion P = u W where W 
is the energy stored per unit length IV^LQA/(L^d) where A is the junction 
area. Thus letting Ug = S/c, one obtains 


a = 





(I(VQ+riw/e) - I(VQ-huj/e)) 


(II. 5. 6) 


If for parameter values one assumes d = lOA, S « 20 , I(Vn+i>uj/e)/A = 

4 2 ^1 

5x10 Amps/cm , I(VQ-fiD/e) = 0, liw/e » 1 Volt, one obtains a c 36 cm" . 

There are two important aspects to this result. Firstly, this a 
must be able to compete and overwhelm the mode absorption for efficient 
conversion. At low temperatures this should be no problem, although Tow 
temperature loss studies should be performed. At room temperature, the 
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mode loss Is too great for this to occur. 

Looking ahead. Fig. (III-8) indicates a mode loss of approximately 
4 5-1 

10 to 10 cm at room temperature so that significantly reduced loss 
must be obtained. For low temperatures (sio®K) three of four orders of 
magnitude reduction is easily obtained so that a can readily dominate. 

The efficiency is then limited by coupling and harmonic conversion. 

Because of the second important aspect even if a is large, one can 
still be left with a low efficiency device. This is clear from Eq. 

(II. 5.1). Since one observes that if I(VQ+tu/e) = -UV^-hw/e), to first 
order I is equal to zero and no conversion to d.c. occurs. This is the 
case for a zero biased symmetric junction device. 

2 

One should thus define the inherent efficiency as the ratio, I 5c'’/" (iiy , 
of the power delivered to load tothe total incident power absorbed. An 
additional mode loss factor of where 3" is the imaginary part of 

the mode propagation coefficient should be included* Thus the overall 
efficiency for the device, excluding coupling is given by 

eV 

c e „, 5 ,, 

(I 

n 

_3 

At room temperature a/(a+23") - 3.6x10 indicating a peak efficiency 
of .36%. In this expression the - I^(Vq) term accounts for the power 
required of the bias circuit. 

It is interesting that n. theoretically at least, can have a value 
of 1. This occurs for a situation for which a 3 " and I(VQ+nhu)/e) - 0 
for n +1 and coupling to the mode is assumed perfect. This in one compact 

statement expresses the essential features for a 100% converter. That 
this is possible is interesting, since, according to this, low frequency 
rectification can be surpassed by preventing multi photon transitions from 
occuring. Furthermore it is interesting '.o note that the above- cannot 
happen for a depletion layer device since phonons inhibit the coherence 

upon which this result depends. 

★ 

While in print it has been observed that recent results obtained by 
R. Drul linger et al , Applied Physics Letter 42, 137 (1983) reports 

an 0 ) ** dependence of the signal to noise ratio of^a point contact diode 

from 1 GHz to 10 THz. It will be shown that and 3 >> a, so that 
these results tend to verify Eq. (1 1. 5. 7). 
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The low frequency small signal limit of Eq. (11.5,7) is obtained 
by taking only the n = ±1,0 terms, taking only the first terms in the 
series expansions for Jq and and approximating differences by deriva- 
tives. Current terms independent of are subtracted since these are 
bias terms. Assuming a » jH'* gives 
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Choosing I7I' s 7(V"^), R"^ ~ I' ^ (300O r .2\l and I - 20 |.A, one 

obtains n s even including the bias power required. These 

results are highly questionnable however, since first of all for 
,2V one needs the large signal analysis which will decrease this projected 
efficiency. Secondly, the results are highly dependent on V^. P.sducing 
this to .IV causes the conversion to beeasily overwhelmed by the bias 
power required. Finally for visible frequencies, for a correct analysis, 
the quantum limit is required since IWe IV, We are presently discussing 
the analysis of this using a more accurate approach for a junction con- 
figuration for which backflow has been diminished. Unless this is done 
it is believed that high conversion is not likely. 

In conclusion, we have carried out a stationary state analysis of 
tunneling and have been able to demonstrate the detailed aspects of 
the interaction of the tunneling with a confined junctiuii mode. The 
final equations for the load current, the inherent transition rate from 
the mode and the expression for the efficiency are simple and provide the 
essential features required to obtain high conversion efficiencies. We 
are continuing our studies of these relationships to first of all investigate 
the fundamental limits of present converters such as p-n junction depletion 
layer devices, and secondly to investigate junctions which will provide 
the highly selective constraints imposed by Eq. (11.5.7) =1. 
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Ill PLASMON MODES IN METAL-BARRIER-METAL STRUCTURES 


For the various tunneling junctions and coupling devices which have 
been proposed to directly convert a free space field to an electrical cur- 
rent through the junction, interfaces between dielectrically different 
media will be important to constrain the spatial extent of the radiation 
field. Such self-sustained electromagnetic interface modes can exist 
provided certain dispersion relationships are satisfied. Since the 
dielectric properties of metals and heavily doped semiconductors at optical 
or infrared frequencies and below are strongly influenced by the collective 
oscillations of free electrons, these excitations are referred to as sur- 
face plasmons. They are T.M. in character and depend upon the presence 
of a negative dielectric coefficient for at least one of the bounding 
materials. 

Some of these modes couple strongly with tunneling electrons due 
to their large dipole moment. For suitable conditions they also couple 
efficiently to an external electromagnetic field. This coupling can be 
accomplished by phase matching using surface roughness, prism couplers, 
gratings, or by confining the excitations to finite structures. The 
various possible coupling schemes are discussed in more detail in Chapter 
IV. In the present chapter, we review some of the basic idealized plasfnon 
modes and then consider aspects of particular importance to potential 
energy conversion devices. In particular, we review some recent numeri- 
cal results which have been obtained. 

III-l The Shape of the Resonance 

An important factor in considering the interaction between electrons 
and surface plasmons is the line-shape of the resonance. This is given 
by the rate of change in energy, as a function of frequency. Consider 
the introduction of a small amount of charge into a surface plasmon 
system. This can be experimentally accomplished by transmitting electrons 
through a thin foil. The total charge p,j which consists of this in ad- 
dition to the polarization charges induced is still zero. For frequencies 
where the electrostatic treatment is valid, the normal component of the 
electric field is equal and opposite on opposite sides of the surface. 
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The latter are related to the charges by: 


onmmx r 

or POOF? QUAUit 


Eq - El = 0 


(III. 1.1a) 


CoEq + EjEj = pf 


Eq ~ P^/(*^g’^6;i ) 


(Ill.l.lb) 


where Eq and Ej are the normal electric fields. The rate of change in 
energy of surface plasmons is given by: 


= Re 
dt 


dV E 


3D* 

3t 


= 1 


dV Re[-iwp| (gq+Cj)*] 
dV Re[-itop^/(eQ+ej)] 


« Im(^) 

£q’*-£1 


(III. 1.2) 


Hence the line-shape is given by Im(---- 7 — ) • For example, for the Drude 

^0 ‘'I 

model dielectric function iQ(l-u)p/c/) , the quantity (cg+f^) has a simple 
zero at w = ui^/VZ. The line shape is thus given by ^(w-mp//?) . A more 
realistic dielectric constant is that for a low-loss plasma. For 
oj » ^ , where C is the mean collision time, this can be written as: 


w, 




^1 = "0 


[(1 - -^) - i 4- ] 


0) 


3„ 

03 t. 


(III. 1.3) 


For frequencies near o3^//?, where Im(- -^;:- -) is significant, the latter 
can be written as: 


"P""' V^I 


Ini(7r4-) “ 




'V^l' (o3-03p/*^)^+(/2/4C)^ 


(III. 1.4) 


that is, a Lorentzian line-shape with a peak at 03 = osp//?, and a full 

width at half maximum of 1/C- In roal metals, an approximately Lorentzian 

cn 

line- shape is obtained from Im( — ~ — ), at least on the low frequency side. 

Cq+Ei 
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On the other side this is usually complicated by interband transitions. 

For a simple metal-air Interface the resonant profile is that of the 
bulk resonance (given by Im(l/Cj)) shifted to a frequency This 

bulk resonance is shown in Fig. (III-l) for illustrative purposes. 
Additional data is available in Ref. (5). Table I shows basic data for 
several of the important bulk resonances. 

1 1 1-2 Idealized Lossless Plasmon Modes for Some Important Structures 
When one discusses the lineshape as in the previous section, the 
dispersion is neglected; that is the dependence of the resonance fre- 
quency on the wavelength of the radiation. This thus gives the resonances 
for the portions of the dispersion curves that are independent of 3. 

To discuss the important aspects of such modes, the full dispersion 
is extremely important particularly since one generally is away from the 
resonance frequency. In this section we consider a variety of ideal modes 
which are basic to tunneling devices and potential coupling devices. To 
deduce these dispersion relationships we initially neglect the influence 
of loss. In the following section this (loss) is considered for two of 
the modes of importance. Recent numerical calculations being carried 
out to more fully investigate the modes of importance within the junction 
devices themselves are also considered. 

The simplest interface is a metal-dielectric plane, as illustrated 
in Fig. (II1-2). The dielectric constant of the metal is c^, while that 
of the adjacent medium is Cg. Consider an electromagnetic wave traveling 
in the z-direction, parallel to the surface. Since a surface mode is 
desired here, the fields are assumed to be exponentially decaying away 
from the boundary, x = 0. With this assumed form, the tangential com- 
ponents of the fields are written as: 

for X > 0: 

Ez = Ag exp(iu)t-ikz-VgX) (III. 2.1) 
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Fig, Ill-l : Optical Properties 
of Aluminum: (a) Reflectivity and 
absorption coefficient data over 
an extended range (Reprinted from 
Ref. 6) 


(b) Frequency dependence of c,, <2 
conductivity o and the energy loss 
function - Im The effect of t 

1.6 eV interband peak on Cj, 
a is obvious. (Reprinted from Ref. 
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for X < 0; 


Aq exp(io)t-ikz+Vj^x) 

itoe, 

_.. _k || _ 1 k r 

wG^ y " vi z 


(III. 2. 2) 


where A-, is the field strength, k is the propagation constant, and 
2 2^ 2 

V = k - 0 ) PE in the respective media. The dispersion relationship be- 
tween (i) and k, obtained by applying the boundary condition of continuity 
of these components, is given by. 



Since Vq and are assumed to be positive quantities, it follows that 
Gq and must be of opposite sign. 

The dispersion can also be written as: 

Of particular interest are those frequencies for which = -Iq, and thus 
k •+ These are the resonance frequencies of the surface plasmons. 

For real metals, g^ is in general complex; k obtained as the solution to 
Eq. (III. 2. 3) or Eq. (III. 2. 4) is also complex, k = P-ia. The decay 
length or range for such surface plasmons, Z = 1/a, which is usually 
very long below resonance [7,8,9] is discussed in the next section. The 
resonance frequencies are still given by cj = -Eq, where is the real 
part of E^. 
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The dispersion obtained above can be expressed in an explicit forni 
using the Drude model for the metal, where the dielectric constant is 
given by; 



where Wp is the bulk plasma frequency: w 
Eq, (III. 2. 4) can then be written as: 


2 . 

P ' Eq"' ■ 


or 
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■peQ = (k 


2 2 
W„-0) 
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(III. 2. 

Equation (I II. 2, 3) and 


( 111 . 2.6 ) 


(III. 2. 7 ) 


The resonance frequency is Wp//? . Such dispersion is representative of 
ideal metals, and is illustrated in . 

For the model above, which indicates the decay into the metal 
and is approximately Wp/ilc^ for w >»- Wp, which for typical metals is of 

the order of (1000 A) , Such thicknesses are not uncommon for metal 
film electrodes in metal -barrier-metal junctions. The interaction be- 
tween the fields of the two interfaces is then significant, and the 
description of the surface plamon mode as a single interface mode becomes 
inaccurate. Both surfaces must be considered in order to correctly de- 
scribe the dispersion. 

The penetration into the air above the metal is similarly given by 
Vq = —(w/peq) for u) < Wp and is thus typically lOX for aluminum through- 
out thi visible region. This however does not take relaxation into account. 

A modification of the simple interface of Fig. IH-2 , the dielectric- 
metal-dielectric structure, is illustrated in fig. III-4 . The fields in 
the dielectric regions are again decaying exponentials, whereas those in 
the metal are a combination of growing and decaying ones. Due to the 
synwiietry of the structure about the mid plane, the fields must either 
be symmetric or antisymmetric. The dispersion for these two branches is 
given respectively by: 
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Fig. III-2. Schematic diagram of an infinite-extent 
planar interface between vacuum and a 
metal. The region above the y-z plane 
is occupied by free space of dielectric 
constant Cq, and that below is occupied 
by a metal of dielectric constant 
Surface plasmon excitations are assumed 
to be propagating along the z-direction. 
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Fig. III-3. Dispersion of surface plasmons on a planar interface 
for an ideal metal, w is normalized with respect to 
ojp, and k with respect to Wp/c. The dash lines 
represent the limits given by w = kc and w - oip/k'^. 
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Fig. 1 1 1-4. Schematic diagram of a vacuum-metal -vacuum 
planar structure, A metal film of thickness 
d occupies the space between x = to x - 
Surface plasmon excitations are assumed to be 
propagating in the 2-dlrection. 
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(^) + tanh(vjd)(-^) » 0 
Vq i Vj 

tanh(v,d)(-~) + (■;-)« 0 

i Vq Vj 


(1II.2, 8) 


For small values of d, the branches behave significantly different from 
the single surface mode. The dispersion for the symmetric branch becomes 
further relaxed, while that for the anti -symmetric branch becomes closer 
to the porpagation mode in free space. The resonant behaviour is 
however very similar in all cases, with the same resonance frequencies, 

For the simple Drude model, the dispersion of both branches are illustrated 
in Figs. (II1-5) and (111-6). 

A closely related but important structure for tunneling junctions is 
a metal-dielectric-metal gap occurring in the junction, as illustrated 
in Fig. (111-7). Here only the antisymmetric branch is of interest, since 
the potential difference between the surfaces along a path perpendicular 
to them is zero for the symmetric branch, and therefore has little effect on 
device operation. From Eq. (111.2.8), the dfiipersion can be written as: 


(^) + tanh(v„d)(i) = 0 (III. 2. 9) 

Vq U Vi 

The field configuration and the derivation of this is given in Appendix (U). 

For very .small gaps the phase velocity is much smaller than the speed of 

n 

light The dispersion is approximately given by: Vq • (Vjt.Q)/{-rjd) . Furthermore, 
Vj, ((3^-u)^yej )*). can be approximated by 3 in some cases or by 


in others resulting in two approximate but useful expression for 3. Assuming 
Vq s 3 since 3 » w/pCq one obtains respectively 


3 = (i“)/d 




(iu.2.l0a) 


,2 . 




;-ej)^d) ; Vj s: r UpV^ 


3 (III. 2. 10b) 


and 


Fig. III-5. Dispersion of the symmetric mode (a) and the anti- 
symmetric mode (b) of surface plasmons on a planar 
ideal metal film. The thickness is given by (wp/c)d = 
w is normalized with respect to w , and k with respect 
to (w/c). The dash line represents that of excitations 
on a single planar surface. 
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Fig. III-7. Schematic diagram of a metal -vacuum-metal planar 
structure. A free space gap of width d occupies 
the space between x = +-^ to x = 4- Surface 
plasmon excitations are assumed to be propagating 


preceding page bank not 


m the z-direction. 
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In the second case, the penetration into the metal is of the same order 
as for the single surface mode. It is llijs latter case which is most 
prevalent as will be seen. These expressions are valid for w Wp, where 
o)p is the plasma frequency. The resonance frequency for the Drude model 
Is again given by Up//2, and the behaviour near resonance is very similar 
to the cases previously discussed. The dispersion using the Drude model 
is illustrated in Figs. (III-8) and (III-9). The circuit properties 
of this mode are discussed later in terms of its important role in the 
coupling of radiation to metal -barrier-metal junctions. 

Although only planar structures have been considered here, the surface 
modes have been considered for other geometries such as the right circular 
cylinder and the ellipsoid of revolution. These could eventually play 
an important role in the coupling of radiation into tunneling devices. 

Recently the plasmon resonance associated with spheroids and ellipsoids 
have become important because of the coupling through surface periodicities 
and roughness. This resonance follows straight-forwardly from the problem 
of a dielectric sphere in a static constant field. If Eq is the field 
amplitude outside of the sphere then the potential inside the sphere 
is given by 

EQrcosG (III. 2. 11) 

For a dielectric constant which is negative and in particular given by 
the Drude model one observes that there is indeed a resonance. This has 
been investigated recently by Liao and others for the enhancement of the 
Raman effect, as well as the coupling to luminescent junctions in periodic 
arrays 

III-3 Propagation Length Limitations Due to Loss 

In all of the idealized structures considered in the previous section 
the neglect of the loss is the most important approximation which has 
been made. For a tunneling converter, the stimulated absorption due to 
the tunneling transition must be much larger than the absorption due to 
the dielectric losses. The latter of course implies conversion to heat 
and hence a decrease in the efficiency of the device. As we see in this 
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Fig. III-8. Dispersion of the anti -syiiimetric mode of surface plasiiions in 
a small planar gap between ideal metals. The width is given 
by (Wp/c)d = 0.10. 0 ) is normalized with respect to top. and 
k with respect to (Wp/c). The dash line represents that of 
excitations on a single planar surface. The phase velocity 
for 0 ) ■> 0 is approximately 0.022c. 





k 


Fig. III-9. Dispersion of the mode for a width of (Wp/c) d = 0.04. 

The phase velocity for cjo 0 is approximately 0.014c. 
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section^ this implies severe constraints on the junction devices, since 
the gap mode in particular propagates a very short distance. 

The most straight-forward means by which to examine the loss is to 
include relaxation effects in the dielectric coefficients of the metals 
and the semiconductors. The simplest expressions by which to estimate 
the resultant propagation lepgths are Eq. (111.2.4) for the single surface 
plasmon and Eq. (III. 2. 10) for the junction mode. The former provides the 
exact propagation distance, however the latter is an approximation, first 
of all, assuming that the same dielectric occurs on both sides of the tun- 
neling junction, and secondly that one is not too close to the plasma 
resonances of the dielectric materials. 

One assumes that the dielectric behaviour is described by 

cjM ' (in, 2. 12) 


where and are the ideal frequency independent dielectric constant 

as 0 ) •+ «>, the metallic plasma resonance frequency, and the relaxation time, 
respectively. One can then obtain approximate expressions for the propagation 
distances for various frequency ranges (a) the low (b) intermediate, (c) 
relaxation and (d) resonance, 
a) Low Frequencies (w£ << 1) 

iM Q 

For low frequencies ej(u)) s iEj^w^C^/w = iOj^^/w where is the d.c. 
conductivity. Then Eq. (III. 2. 4) gives for the propagation length of the 
single surface mode k" (k=k'-ik") 


k •' = ik^ (a„y^ )-' (III.2.13) 

and k' - k = . Thus k"/k' = i(k)(o,^, /^)’^ a long distance 

/ ‘" o 

(Ojj^ = 10^, thus (k")”^ s (lO^)X^ cm). The penetration depth into the 
metal reduces to the skin depth and the penetration into the air is given 
by Yj = ^ ^ rad/sec. 

The metal-barrier-metal or gap mode propagation distance follows from 
Eq. (III. 2. 10b) and in the low frequency conductivity limit gives for 
k = k'-ik" 
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(III, 2. 14) 


where the .38 and .92 are respectively cos(~) and sin(|-). The most 
striking characteristic is the uF dependence for both k* and k". 

Thus the real part of k in this region is approximately 2.4 times 
the imaginery part and the ratio of k'/kQ * (2.42)k"/kQ « l/(d^wpQO|^^)^ 
s (6/d)^. The propagation distance is thus of the order of which 

is about 10“^X at 10^^ radians per second, an extremely short distance 
(sl8pm) 

b) Intermediate Frequencies 

At higher frequencies for which the real part of the dielectric 
begins to contribute significantly the expression (2.14) for k of the 
gap mode is replaced by 

sin(J-wjy4)) (III. 2. 15) 


k 


where « 1 

21 This shows that the real part continues to increase as approximately 

but that the imaginery part increases less rapidly with increasing 
frequency. is long enough, k" can in fact decrease and thus the 

propagation distance can remain constant or decrease slightly. For room 
temperature this is not expected to be observable. For the single 
surface mode expression (2.13) given above is still valid in this fre- 
quency region, 
c) Re'axation Region 

This corresponds to frequencies for which the real part of the dielectric 
becomes comparable to the low frequency conductivity contribution or is 
the dominant contribution. This is expected to be the most important 
region with regards to convertors since this, aside from interband 
absorption, is applicable for the optical portion of the spectrum. The 
dielectric coefficient can be approximately written as 
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(III. 2. 16) 


Gj(w) 


- 0 ). 


+e. 


m 


0)^ 


w 


+ i 


m 


3r 

0 ) h 


ni 


where it is assumed that » 1. 

Then one obtains from Eq. (III. 2, 4) for the single surface mode 




nr ■'m' 


(III. 2. 17) 


where c is the speed of light, Consequently the propagation distance is 
roughly a constant of the order of 3 pm for a relaxation time (r ,) of 
the order of 5x10 sec. Since for metals Cq = propagation distance 

is roughly given by the relaxation time in wave numbers, 

In this region of the spectrum, the penetration of the mode into the 

dielectric is given by Vn ' ’Jt" * frequency limit, and 

^m m ^ 

the penetration into the metal is given by Yj ^ 

For the gap mode Eq. (111.2. 10) with this approximate value of 
Ej(m) gives 

^ (111.2.18) 


One expects for » 1, the real part k' to be proportional to w 
and k" to be approximately constant and given by 


(k') = 

~ 1 1 1 ~ 

”(m E^ph^ d'* 2(Gf;„) " 
' m ni 0' ' "m' 


1 

( 10 ’^) 



(III. 2. .19) 


assuming k^ = (w^^/e^Pq) s .lo'^ cnf^, d lO'^ cm and t 10"^^ sec. To 
observe the constant k", Cj^ must be long enough to allow this region to 
be observable without having the wing of the resonance peak dominate the 
behaviour. 

d) Resonance Behaviour 

For frequencies for which the influence of the resonance comes -into 
play, one can assume in Eq. (III. 2. 4) that in the numerator is given 
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predominantly by the real part and that this is approximately given by 
-Gq, One then obtains 


3 2 

2 . “r^O^O 


1 




r' 1 

(a)-w>l7rx-n4“) 

“ni 


(111.2.20) 


which displays the resonance profile shown In Eq. (1,4) where = 




e_ i 
m 


c„+e, 


) . This expression can be used to obtain the propagation 


m’"0 y C i 

distance. At resonance this is oj (— ^ -- Q - z - -■ -- ) . Eqs. (2.10a) and (2.10b^ 


" VS 


although not expected to be strictly valid, show that the gap mode displays 
a simillar resonance determined approximately by Re Cj(w) = 0, that is, 
close to u)j^. These two expressions (2.10a and 2.10b) become 


3 


1 ^0 2 

i - ^ 
m 




( in , 2 . 21 ) 


111 


"m 


and 


1 ..2 *'0 


ii) 


} ^0 ^ 


1 


m 


/ 2. 2 ,/Svs 


(III. 2. 22) 


1 

The first predicts that at resonance 3' = 0 and 3" = ’p(S<‘',n) and 

Q (iJ j I 

that a Lorentzian profile is obtained. The second gives 




-j(-i)^ So that 1 3" I s 3' and above the resonance the character of the 
m 

mode changes. It is this latter behaviour which is found by solving 
numerically. 

In Fig. (III-IO), we have shown propagation distances for the open 
surface mode and the junction mode using numerical evaluation of expressions 
(2.4) and (2.10b) for the case of Indium (hw^^ = lleV, = Cq). One 
observes that the propagation distance of the open surface mode is tolerable 
throughout the visible portion of the spectrum at room temperature (" 50. )iiii) 
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The gentle rise throughout the visible is due to the resonance which is 
more pronounced for the single surface mode. 

Below resonance the relaxation results are extracted for w 
Eq. (2.22) gives Eq. (2.17) as is seen by expanding the denominator. 

Eq. (2.21) gives 

1 Eq j^2 j ORIGINAL PAGE S9 

^ ® B F" “T*' OF POOR QUALITY (I II. 2. 23) 
m ■'( 1)111 

which provides a quadratic dependence on w as might be expected in the 
blue portion of the visible spectrum. 

For the case of a metal there is considerable advantage to cooling 
since the relaxation time increases with a decrease in temperature, and 
this is proportional to the propagation distance. In Fig. (IlI-lO), we 
list the propagation distance for the open surface and junction modes as 
a function of the relaxation times and show that the propagation distance 
is enhanced considerably. There may be additional inherent advantage to 
running the junctions cooled to obtain such relaxation times, One would 
expect the tunneling characteristics to be more ideal as the temperature 
is lowered and thus It would be expected that the transition probability 
for the absorption processes should increase slightly with a drop in 
temperature. 

There does not appear to be much data available in the literature 
relating to the temperature dependence of the relaxation times of metals. 
However the data which is available does indicate that the relaxation 
is in fact phonon dominated and can thus be considerably lengthened by 
operating at reduced temperatures, and in particular the photon pump 
processes might be able to exceed thermal i zing processes reasonably well 
but that the propagation distance of the junction mode is extremely 
short (s .5pm). More work needs to be done to investigate such relaxation 
processes at lowered temperatures. 

III-4 Numerical Solutions 

We have in addition to the above analytical approximations to the 
junction mode, numerically obtained the dispersion characteristics of 
the gap mode. This is appropriate when the dielectric material on one 
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side of the tunneling barrier differs from the material on the other 
side of the barrier or for more complicated barrier structures that will 
probably .leed to be considered in the future. For the case which we have 
explicitly considered, the dispersion expression is given by 




fl 







“ci 


fnl 

-J 


>- 

ill 


-Yid 


(III. 4.1) 


where the transverse propagation constants Yg* Yp Y„i are given by 
Y = [3^ - ^ ® S’ *^b’ '^s dielectric constants 

of the respective regions shown in Fig. (III-7). They are taken to be 
e^(w) = e^(l-a3^/((D^+iw/r,^,y), Cj = and Gg(cu) « ^‘^(l-Ws/Cw^+iw/r,^) ) , 

where to , ^ , w_ are the relaxation times and plasma frequencies of the 

d Ml PM 

semiconductor and metal respectively, and Cj, are the ideal frequency 

independent dielectric constants as w in this model. Fig. (111-11) 

illustrates an example calculation for a dispersion curve using this 

numerical procedure. For this particular case was arbitrarily chosen 

for Aluminum (tiw_ s; lleV) and the behavior is dominated by to. and i;.. 

A relaxation time ® 4.0 x 10"'^'’ sec. was chosen as v^ell as hto. = 

20 ^ -11 

.35eV (heavily doped semiconductor (- 5x10 cm )) and ^ = 1 x 10 sec. 

2 3 “ 

The relaxation times are approximately 10 to 10 longer than expected 
but should be achievable by cooling to reduce phonon domination. These 
results should thus provide an indication of typical behavior expected 
of cooled metal -barrier-heavily doped semiconductor devices. One observes 
a much reduced loss. Fig. (III-12) and (III-13) show siiniliar calculations 
for Al-AlgOg-Ag and AI-AI2O3-AI junctions respectively. Fig. (III-12) 

also indicates the various regions of the dispersion curve. 

The log-log scale clearly illustrates that at low frequencies one 

observes the dependence as suggested by Eq. (III. 2. 14) with the real 
part being greater than the iniaginery part by approximately a factor of 
2.5. At higher frequencies for which the real parts of the dielectric 
coefficients contribute significantly along with the imaginary parts 
(relaxation region) the dispersion deviates from such a dependence. One 
observes an w dependence for 3' and an almost frequency independent 3" 

(as expected in the relaxation region). 
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One observes from Fig. (III-12) that for hw ~ leV, one is in this 
intermediate range between the relaxation behavior and the resonant 
behavior, A value of 0“ of approximately 2x10^ or a propagation 

distance ^ .0 ym for the c 6x10"^^ sec can be expected. This is a 
reasonable propagation distance with respec; U the absorption which could 
be anticipated from relatively nonlinear device; (see £q. II. 5. 6). This 
propagation length, also in good agreement with Eq. (2.18), could also 
be lengthened by cooling. Since the relaxation time decreases linearly 
down to the Debye Temperature, jej, and as T^ thereafter Eqs. (I II. 2. 18) 
and (111.2.22) indicate that by cooling to approximately 10° K the pro- 
pagation distance is increased to at least 200 ym and thus the tunneling 
transition rate can easily compete with the plasiiion propagational loss. 

There are two other parameters which are significant when one is 
concerned with maximizing the transition rate associated with the tun- 
neling phenomena. The first is the confinement factor or the fraction 
of power that is flowing within the tunneling layer. This can be large 
and can exceed unity since the negative dielectrics cause reverse power 
flow. For the region of interest. Fig. (III-ll) indicates that this con- 
finement factor C is close to unity. The power confinement is excellent 
although the energy is principally stored outside of the tunneling 
barrier. This is of course due to both the large negative dielectric 
coefficients as well as the longitudinally directed electric field com- 
ponent of the plasmon mode. The latter becomes important near the 
resonant frequency of the mode. The question of power and energy confine- 
ment needs more consideration. 

The second parameter of importance is the slowing factor, the ratio 
of the speed of light to the energy velocity associated with the mode 
and defined in the first chapter. One observes from Fig. (III-ll) that 
this exhibits divergent behaviour at low frequencies due to the influence 
of the conductivity but that for the frequency range of interest it is 
relatively constant with a value of 200 confirming the number used in 
Chapter II to estimate the absorption coefficient due to the converting 
transitions. We are presently in the process of calculating the dispersion 
characteristics utilizing more accurate values for the dielectric para- 
meters as a function of temperature. In addition, in the visible region 
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of the spectrum, for most metals since interband transitions add at least 
a small amount of additional absorption, these are being included in the 
calculations. Finally more complicated multilayer tunneling structures 
having the form of metal -semiconductor-tunneling barrier-semiconductor- 
metal might have to be considered and present more difficulty in calculations. 
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IV COUPLING BETWEEN FREE SPACE RADIATION MODES AND JUNCTION SEW MODES 
IV-1 Introduction and Overview 

In order to extract energy from free space radiation modes, such 
modes must be coupled with low loss to the well confined surface electro- 
magnetic junction modes. The coupling devices required to accomplish 
this must basically guide the incident field from dimensions of the order 
of microns to dimension of the order of Angstroms, a four order of 
magnitude change in the confinement distance. It is quite likely that 
this coupling will involve the intermediary of a single surface plasmon 
wave. Thus the coupling will involve two steps; the coupling of free 
space modes to single surface SEW's and subsequently the coupling of 
these to the highly confined junction modes. 

Thus one must basically consider three problems involved with the 
coupling. 1) The coupling of the incident field to a guided single 
surface electromagnetic mode; 2) the coupling of the single surface 
electromagnetic mode to the confined junction mode, and 3) the coupling 
of the junctionmode to the electron tunneling current. 

In the present chapter we wish to consider the first two coupling 
problems and report on the progress which has been made to alleviate 
them involving both our own efforts and those of other groups. 

An overview of various coupling schemes which have been reported 
is shown in Fig. (IV-1). In the first three the technique is to obtain 
an overlap of the free space radiation modes propagation dispersion relation 
ship and the dispersion relation for the plasmon modes. In the last 
Fig. (IV-l(d)) the diffraction pattern of the plasmon mode from an abrupt 
termination is utilized as a radiative aperture for the coupling. The 
first is analogous to integrated optical coupling (the Otto configuration^^^ 
The prism refractive index n is large enough to allow the radiative 
mode phase velocity ^ to exceed the plasmon mode phase velocity. A finite 
angle 9 is thus obtained for an intersection of the two dispersion curves 
and hence the possibility of coupling. In the second technique (the 
Kretschman configuration^^^ ) , the plasmon mode on the bottom surface 
of the metallic film is evanescently coupled to the radiative mode 
utilizing the exponential decay in the metallic film. The length 
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confinement coupling of Fig. (IV-l{c)) is easily understood in terms of 
the uncertainly relationship AkL > 2ti. Thus Ak ^ (2 ti/L). This smears 
out the dispersion curve and causes an overlap with the light line and 
hence coupling. The lobe pattern results from the detailed spatial 
frequency Spectrum of the confined mode, which has a total width of Ak. 

We shall consider this in detail subsequently. One observes that if the 
length is infinite the SPW does not radiate in this configuration. Various 
combinations of the above schemes are of course possible, in particular 
length confinement and evanescent wave coupling. Length confinement 
results in side lobes on the evanescent wave coupling characteristics. 
Similarly, the boundary conditions associated with the end coupling 
characteristics will influence the length confinement mode pattern. The 
prism coupling schemes to date have not achieved adequate coupling for 
efficient energy conversion (on the order of a few percent)^^^, although 
it is felt that optimization could increase the coupling to acceptable 
levels. 

Prism coupling using end radiation has achieved the highest coupling 
efficiency of any method' \ Dsses are due to parasitic scattering and 
inefficient mode pattern matching. It is felt that this method is 
potentially promising for efficient coupling, as optimization reduces 
to edge preparation and mode matching by a suitable lens. It is also 
only slightly frequency dependent. 

Length confinement as an extension of antenna theory ^^’^^takes into 
account the high frequency dielectric behaviour. It thus should prove 
possible to have well designed antenna arrays for matching to the 

radiation field and thus these should be potentially efficient couplers. 

(7l 

Grating coupling'^' has also been used for coupling single frequency 
radiation to the SPW, but suffers from low efficiency due to reradiation 
and coupling to higher order modes. Grating fabrication, being a planar 
process, is however directly compatible with thim film technology which is 
advantageous for the construction of thin film energy converters and 
progress has recently been made on increasing the efficiency at submilli- 
meter wavelengths^^^ . Gratings are producible through direct photolitho- 
graphy, laser interference, or laser photodeposition^®^ . The latter method 
shows promise as the SPW coupling is optimized, i.e. the grating profile 
grows in such a manner as to maximize the SPW intensity in the presence of 
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a radiation field. It has been shown that in the presence of dielectric 

{Q\ 

losSi total absorption of the monochromatic incident wave is possible' ' 

with a suitable grating profile. As in the evanescent coupling schemes, 

however most of the energy is absorbed in the materials, hence actual 

mode conversion efficiency is low. 

The strong dependence of coupling angle on frequency makes periodic 

gratings unsuitable for wideband couplers. Related methods, e.g. aperiodic 

gratings and surface roughness^^^\ show some promise in this application, 

but again coupling efficiency is fundamentally low. Nonlinear optical 

surface studies^^^^ and investigations of laser induced surface gratings^^^^ 

have recently been reported and should result in development in this area, 

In addition enhanced plasmon mode coupling to an array of junctions at 

(13l 

visible wavelength has recently been reported' . 

In the present chapter of this report we consider primarily length 
confinement coupling since this is the area in which our concentration 
has been. We briefly consider taper coupling which can provide efficient 
conversion from the single surface plasmon mode to the well -confined gap 
mode and which we are continuing to pursue. In the final section we report 
on recent theoretical investigations involving grating coupling to the 
single surface mode. 

IV-2 Length Confinement Coupling on a Planar Structure 

Although a true surface mode on an infinite surface does not inter- 
act with the radiation field because of a lack of phase-matching, for 
a finite structure such a coupling is possible. In this section we show 
that such surface plasmon excitations on a linear structure properly 
reduce to conventional antenna characteristics at low frequencies. It 
is also demonstrated that such plasmon excitation provides efficient 
coupling between the radiation field and a metal tunneling structure up 
to a frequency at least as high as the plasma resonance. We have calculated 
the radiative patterns of such length confinement coupling characteristics . 

The ideal structure of interest is that shown in Fig. :(IV-2(a)). For 
simplicity we consider a planar geometry for the plasmon mode on the 
boundary. Such a structure supports surface modes in regions A which are 
coupled to the radiation field, and an antisymmetric "gap" mode in region 
B which can be efficiently excited by conversion from the former mode. 
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The basic surface plasma mode is that propagating along an Infinite 
plane surface separating two regions with dielectric constants Cj and 
Cg respectively, as shown in Fig. (IV-2(b)), From Chapter III, the 
dispersion of such a mode at frequency w is given by: 

(Cj/Vj) + (C 2 /V 2 ) ”0 (IV-1) 

2 2 ? 

where v in each respective region is defined by v - k - w pc, and k 
is the propagation constant for the mode.* For real metals c is in 
general complex due to losses, and k is complex: k = 3 - la. The range 
given by l/2a is usually much larger than a wavelength, except near resonance 
as we have seen. 

In the region of the gap (region B), which is formed by two plane 
surfaces, of the two non-radiative modes which can exist the anti -symmetric 
is of interest. It has been shown that the dispersion for infinite 
extent of the surfaces shown in Fig. (IV-2(c)) is given by: 

(Cj/vj) sinh(v|d) t (e^/vl) cosh(v^d) = 0 (IV-2) 

as shown in Chapter 3. The corresponding z-polarized component of the 

i(wt-kQx) 

electric field is given by E^Cx.z.t) = Eg(z)e ^ /2 + c.c,, where 

k is the propagation constant for the gap mode obtained from Eq. (IV-2) 

9 2 2 2 

(as above, v = kg - to pc), and the complex amplitude Eg(z) is 

Eg(z) = E^ cosh(v|z)/cosh(v|d/2) , * z ■' (IV-3) 

The magnetic field component associated with this surface electromagnetic 

itoCpE 

mode is y-polarized with amplitude (from — ^'■) H^(z) is given by 

y iKg g 

(OCoE . A ri 

Hg(z) = ~ cosh(v|z)/cosh(v|d/2), "2 < z < '♦f 

= cosh(v2z)/cosh(v2d/2) (IV-4) 

Although in general kg cannot be expressed in closed form, for 
frequencies such that the penetration depth into the metal is much larger 
than the gap width, which is generally true for M-B-M structures. 


* Note; in accordance with Fig. (IV-2(c)) we use l '2 = Cq and Yq " Y 2 
A superscript g is used for the gap mode. 
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v| a kg, and vj « /-ywcj and one thus obtains the approximate results 
obtained earlier and shown numerically to give at least quantitatively 
reasonable results: 

kg = (k^ - (IV-5) 


This corresponds to the relatively flat region for 3"(u) from Fig. (1 11-11). 

The characteristic impedance of the anti -symmetric mode for a width 
2b in the y-di recti on is given by the ratio of the voltage across the gap 
to the current induced on the upper electrode; that is 


Zg - 2/ Eg(z)d2/(mC2Ei/kg) 

tanh(v|rj/2) 


(IV-6) 


For the length confinement coupling of this section one would expect 
continuity of the direction of the y-component of the magnetic field 
across the gap as indicated on the upper part of Fig. (IV-2(a)). Since 
the x-directed component of the electric field in region A becomes the 
z-component in region B (the gap) this component is in opposite directions 
at z = and z = ■+-|. The z-directed component of the electric field in 
region A can be calculated from these two components since (io3C2E2 " 

= -Y 2 Hy) where Y 2 = ^ surface propagation 

vector. This camponewt one would expect to be continuous across the gap. 

With rffpect to relatively long wave-length surface plasmon wave in 
regions A there is a step' discontinuity in the x-coniponent of the electric 
field due to the anti -symmetry. Thus one can express the field as 


E^(z) = Eq sgn{z)(e"‘'*^l^i) 


Similarly 




where it has been assumed by symmetry that 
E^(z'O) a 0 


(IV-7) 


(IV-8) 


(lV-9) 


* A more general expression for the characteristic impedance valid for 
all d is given by 
k 

Zg = - (tanh(vtd/2) t (v^d/2) {cosh^{v 2 d/ 2 ) ) ) 

^ b 


Similarly one can have plasma waves propagating in the opposite 
direction. The relative amplitudes of the two are determined by the 
boundary conditions at the ends as well as the power radiated. The 
boundary condition in general results in an end coupling condition which 
can provide strong radiation. We have not studied this aspect of the 
problem, but plan to do so in the future since it appears promising and 
a1 so is more broadband than length confinement coupling considered here. 

The latter ideally does not allow end radiation since one generally 
assumes an ideally terminated open or short circuit, All of the power 
propagating toward the end is thus reflected back. This is accomplished 
by terminating the antenna in an abrupt change in metal thickness (Fig, 
IV-3). The impedance changei even in the visible, is large enough to 
establish nearly 100?^ reflection. is forced to become very small 
far away from the plasma resonance for b) the short circuit termination, 

For a) the "open circuit termination" E 2 (z“.'L) 0, although this can 

provide strong coupling from th(?i end of the metallic antenna. 

Using zero end conditions for E^(z) gives for the fields outside 

of the metal 

E/z) = Eq e (IV-lOa) 


Hy(z) = ^ (IV-lOb) 

E^(z) = ^ e (IV-lOc) 

The fields inside are the same with Eq replaced by in E^(z) and Yq replaced 
hy -Yj to guarantee continuity of and E^(z) and pfovide a jump in E^. 

By virtue of the termination one expects the field pattern given by 
Eqa (T0a,b,(:) to radiate. To deduce this radiation pattern one can obtain vector 
potential A far from the antenna by replacing the surface mode with a complex curr« 
The current is given by ia}(c-rQ)E for the fields shown. Since one expects 
Ej,(z) to be the influential term, we consider only this component. Thus 
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lV-3 Ideal antenna terminations, a) Open circuit b) Short circuit 
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X < 0 


(IV>U) 






for the respective current densities in the two media. 

The radiation pattern is obtained by solving for the vector potential 
in the Lorentz guage so that only currents determine A and hence the 
magnetic field. One observes that the above relationship gives a long 
wire antenna pattern at low frequencies. Assuming a width Ly and that 
the decay lengths Yq Yj are much less than a wavelength, one can 
obtain an ideal terminated line current source 


+00 


u)L 


“ / J7(z)dxL^ = V’ c 


•oo 


^ (c 


o)EoK^' 


•ik Iz 


-e 


•ik(2L- 


] Amps. 


= Ine'^*^*"i sin k(L-lzj) 


(IV-1?) 


where Iq has been defined by 


2o)L. 


^0 


(IV-13) 


^0 ^ "k Yj 

In the low frequency limit this reduces properly since = -j 

"fQ £qU) 

Vj s /37(S , k s o)/c and ” s the single surface dispersion 

1 


-£L 

w' 


relations. One then finds Iq s -2Lya6 (Amplitude of the z-component 

Yq 

of the field (jj^ Eq)), which one would expect. Here 6 is the classical skin depth 

From the vector potential one can calculate both the radiation fields 

and the total radiated power. The electric field in the 0 direction is 

given by 1 ^ 

, AI 'o -ikL , 2„(cos(kL)-co5(koLcose))(-^) 

En(e.r,(|)) = /^ — e sin 0 — r — 


O-(-f)Wo) 


(IV- 14 ) 


and the radiated power per unit solid angle by 
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(-2qL) 


— — 5r[cos(k^L cos 0) - cos(kL)] 
l-(kQ/k) cos'^e ^ 


(IV-15) 


The two extreme values of (kQ/k) provide two limiting cases for the 
coupling characteristics. For low frequencies (kQ/k) s l and ai « 1. 
The radiative pattern reduces to: 


dP 


R\ _ n 




Stt 


2 *0 


cos(kQl cos 0) - cos(kQL) 


sin 0 


(IV-16) 


This is as expected is identical to that Of a standing wave long wire 
antenna, and is shown in Fig, IV-4. The more relevant case for infrared 
and optical wavelengths at room temperature is the high frequency limit 
for which (kg/k) < 1 and aL 1. Here the radiation pattern as illustrated 
in Fig. IV-4, is given by: 


A) =, -Tl. ,2 


3271 ‘ 


•^0 

2 

sin 0 

k 


l-(kQ/k)^cos^o 



(IV-17) 


The most apparent difference from the low frequency pattern is the lack 
of lobe structure due to the short propagation length of plasmons at 
high frequencies. This aU-o removes the dependence on L of the radiation 
pattern, If 1/3 < L < l/a, lobe structure is expected to be observed, 
but this is in general different from that given in Eq, (IV-lfi) since 


3 > kg. 

The coupling of the plasiiion surface mode to a laser cavity can be 

calculated by a modal expansion using the laser cavity modes. To do 

this consider the electric field of the Gaussian mode in the far field 
TiuiSn 

(Irl Zq = . The far field of the lowest order mode, Eqq is given by 


1 -(x^+y^) 

^00 

where the normalization is such that 


-jk • r + jwt] ( IV-18) 

2 

|Eqq 1 dx'dy' = 1 and w^(z) = 
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(b) high frequencies 
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Radiation patterns at low and high frequencies for the total 
radiated power* 
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2 

and Zq = TTO)Qn/A. Similar expression are valid for the higher order modes 
^mn ' 

We approach the problem from the point of view of a transmitting 
antenna exciting the laser cavity as illustrated in Fig. (IV-5). One thus 
expands the electric field excited in the laser cavity by the antenna 
in terms of these orthonormal modes of the laser cavity. This is of the 
form 


E - I 

where we are particularly interested in a^Q, the coupling to the lowest 

order mode. As long as this is done in the far-field of both the laser 

cavity modes and the antenna radiation field, the coefficients a „ are 

mn 

easily deduced. This is assumed to be the x'-y' plane of Fig.(lV-5). The 
laser field is assumed polarized in the plane of the antenna. The 
coupling coefficient of interest Bqq is then given by 

+00 

; dx'dy'E5(,Ej|(x’.y') (IV-20) 

—00 


where Eg(x',y') is the electric field component of the radiation field of 
the antenna. 

This integral can be easily carried out when the area of the laser 
is much smaller than the lobe size of the antenna and assuming that the 
coupling lens is such that the confocal spot is at the antenna as shown. 
Then r H R, e s |0| and (}) s Q and one can separate agg into a product 
of an integral over the laser field and the field of the antenna. Thus 


‘00 


+00 

= [/ 


dx*dy* 

,w{z)a)(z) 


(2Tir 




z) e' 




]E„(R, 


{ 27 i )‘ 


u.(z) 


In the far field w(z) = 


and ojg is the beam radius at 


Eg(R.]0l,4.) (IV-21) 

^0^ RA 

where n is the refractive index 

Zg 7TWgn 

the antenna. Thus 
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IV- 5 Coupling of a Gaussian beam to a long wire antenna lobe. The 
antenna is oriented along the z-axis at R=0 and coincides with 
the focal region of the Gaussian beam. 


ommm. 

= Ofi TOOK QUAim 

The power coupling to the lowest order mode of the laser cavity is 


thus 


^c ' ^l®00 




2 dP 


^c " ^dQ^^ 

is given by Eq. (IV-14) and jb 


where 


dn 


(IV^23) 


8c' 


00 ' 


2 2 2 

(OqW n 


As an illustration of this, the low frequency short dipole of length 

Lh«(^) = jLL2 


(Attc)' 


n / *jj where Iq is the current amplitude. 


Thus 


uwo 


(IV-24) 

Since the radiation resistance of such a short dipole is R = / ^ (~)^ 
for a conjugate match of the antenna the power from the generator is 
(IqR/ 2) = P^. Thus from the above 


Pc = ( 


Ssin^e 


■) p, 


IV-25) 


where S is the confocal area of the laser at the antenna ttojq. By 

integrating over all angles (f sin 0 dG = 4) one obtains the usual result 

b 

for the effective area for the short dipole antenna. 

To complete the calculation of the coupling of the gap mode to the 
laser field the value of Eq must be expressed in terms of the gap field 
E.j by considering the coupling between the open surface mode and the gap 
mode. Once this relationship is established one can apply the principle 
of reciprocity which demands that the coupling be equal for receiving 
and transmitting. Ideally one would desire one hundred percent coupling 
of the power from the open surface to the gap. One has Power = Re/jE^H*dx 
for the gap mode propagation in the x-direction. For the components E = 
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H for the single surface mode and E_ ® for the gap mode, thus 

u)£ 2 y ^ y 

the one hundred percent coupling constraint implies 

k / III 1 2 


ok 2 

/ m l^clx = -S- / iH^t d2 
a)e« * I y' 0)c« ^ ' V* 

‘^Single 

Surface 

we. 


'Gap 


(IV-26) 


WCp 

One also has the relations H,. = E- for the gap (Eq. IV-4) and 

y '^q 1 

^2. _ Y2 




Ea = 


iwe. 


H for the single surface mode. 


Since the power flow is primarily outside the negative dielectric 
regions, these integrals give for the field amplitude constraints 


WEo o 

2WU 1^0^ 


[g w 2y2^^2 w ,.g 

Thus in the gap c-5n be expressed in terms of Eq. 

Such a coupling from the gap to the antenna requires a taper design 
which can in principle be accomplished for any antenna length. We are 

continuing to investigate this possibility theoretically and experimentally. 

If it is assumed that Kirchoffs Laws are valid for the termination 
(no significant radiation at the junction antenna boundary), then the 
junction coupling problem reduces to a transmission line analysis 
(Fig. (IV-o)). Letting Z (0) be the impedance of the single surface mode 


Q wCod o k_d 

|H„I =-xr |E,-| ' iHj 


(IV-27) 


at 2 = 0, then the relationship between a current source I^e 
at the gap and is 


jwt 


^s ' '^^0® 


-jkL 


sin kL(Hy^ 2^(0)) 


appl ied 
(IV-28) 


where is the impedance (assumed lumped) of the gap. Thus Iq can be 
expressed in terms of the driving current amplitude I associated with the 
gap through this relationship. 

The impedance of the gap Zg(0) is given by Zg/(i tan(kpt)) which for 
a short gap is Z /ik £,. Thus using Eq. (IV-6) one obtains 


1 


1 


be 


1WC/ 


iwC. 


where C is the gap capacitance. Z (0) the impedance of the antenna is 
9 5 

similarly given by 


Zs<0) = ± 


1 

s 


itan(ke) 


(IV-29) 
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IV-6 Illustration of reciprocity for evaluating the receiving charac 
teristics of a laser illuminated antenna. 


- 87 - 




Using Eqs. (14.22,28) one can thus write 
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(IV-30) 


where Z the coupling constant is given by 

JkL 


/H ... 4 Iff.n 

» s' ' \ u 


jsin(kL)(l-fjr^) 


(IV-31) 


P coskL-cos(kQLcose) 

where f(6,L) is sin‘^6( ^ — 5 — ) for the case of transmission 

l-(kQ/k)‘^cos‘^e 

line coupling to the gap mode. As a consequence of reciprocity, one thus has 


V, = (Z)/|a„o (IV-32) 

as the voltage induced across the gap by the laser field, 
since /| ,00 is the magnitude of the magnetic field within the laser 

cavity and also because of normalization corresponds to current. Knowing 
the value of the voltage-current characteristic of the junction allows 

one to thus calculate the rectified signal induced due to nonlinearity. 

The above basic approach shows that antenna theory is easily extended 

to higher frequencies. These results form the basis for our numerical 
estimates of the signals available from the diagnostic heterodyne experiment 
considered in a latter section. In the following section we also verify 
that for the proper taper one can expect nearly one hundred percent coupling 
between the open surface mode and the gap-mode. 

Before leaving this section we briefly consider the gap mode in the 
presence of electron tunneling and some initial calculations of induced 
optical voltage using a more conventional antenna approach. For 
typical junction areas and distributed tunneling conductance Gj(n /m ) 
the gap mode is only slightly perturbed by the tunneling. The potential 
amplitude induced across the gap by the incident field can be obtained 
byi considering power flow through the junction for which the relative 
phases of the waveforms need not be considered. For the ideal structure 
of Fig. IV-la assuming for simplicity symmetrical excitation at the 
edges (x=±a) a standing wave gap mode is excited. The distributed 
tunneling conductance can then be replaced by an equivalent lumped 


a -2k!!(a-|x|) 

conductance = 2b/ Gje ^ dx) in parallel with the impedance 

-a 

of the gap mode as seen from / ±a. By symmetry the gap mode current 

must have a null at the center (x=0). The equivalent impedance of the 
gap observed at the edge is then Z » Z^/i tan(kga). This reduces to that 
of a lumped capacitance of the same dimensions at low frequencies since 
tan(kga) = k^a and Z 'v l/iw(2cQab/d) (using Eqs. (VI-4) and (VI-5)), but 
approaches the characteristic impedance of the anti -syimie trie mode at 
higher frequencies (tan(kga) -i). 

The equivalent circuit consisting of a voltage source with an internal 
impedance driving the above combined gap mode impedance has been used to 
calculate the optical voltage excited by an external radiation field. 

Using conventional antenna theory, the voltage source driving the antenna 
has an amplitude where is the incident field amplitude aQg/C^^irwQ) 

pjkL 

and effective antenna length which for a 

short antenna in the maximum gain direction is given by the shorter of 
the actual length L or the plasmon propagation length 1/a. The internal 
impedance given by the series combination of the radiative, dissipative, 
and reactive impedances is included in 2^(0). The results of such a cal- 
culation carried out by Siu for the optical voltage across an effective 
tunneling resistance (1/Ggff) of 500 ohms in both AJt and W structures is 
illustrated in Fig. IV- 7. For these computations the experimental frequency 
dependence of the dielectric coefficients has been included. Since the 
dissipative and reactive power in the metal become significant, the optical 
voltage decreases as the frequency increases in the region where the 
plasmon propagation length 1/a is much longer than the actual length. 

As the plasmon propagation decreases below the actual length, the 
dissipative and reactive impedances decrease, and since the effective 
antenna length also decreases accordingly, the optical voltage is approxi- 
mately constant. An interesting phenomenon is observed at those frequencies 
where the interband transitions are strongest. Dus to better impedance 
matching from the increased reflectivity at the gap, the optical voltage 
at the surface resonance in AX,, since propagation of surface plasmons 
is inhibited. At frequencies above the plasma resonance these structures 
behave as dielectric antennas, and the optical voltage developed can be 
significant. 
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IV-7 Theoretically predicted frequency dependence of the optically 
Induced voltage. 
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The decrease in optical voltage as the fr^queiicy iincreases frow the t 

infrared to the visible is much less dramatic than that predicted by conven- ^ 

tional antenna theory. The ratio of the optical voltages at 10.6 sa and 

o 

at 6328 A is approximately 8, both for At and W, whereas the ratio! of 
the respective wavelengths is approximately 16. It Is found 
that in the relaxation region the reactive power in the metal becomes 
more important than the dissipative pow'.T. This indicates that by 
properly matching or tuning, the optical voltage can be substantially 
enhanced. 

More extensive investigations of the behaviour of the optical 
voltage or power coupled into the junction at visible frequencies are 

I 

presently being pursued on the basis of Eqs. (lV-31) and (IV-32) to 
calculate directly in terms of the laser field. The important impli- 
cations of an extension of antenna theory through the use of surface 
electromagnetic waves are, however, apparent. 

IV-3 Modes Along a Corrugated Interface - Grating Coupling Basic Relations 

To couple energy to and from a surface plasmon mode at a metal -dielectric 
interface, one might consider the use of a corrugated interface region. 

Coupling between the bound mode and free space modes is achieved through 
Bragg scattering from the periodic structure. In this section, the 
equations describing the field amplitudes along a sinusoidally corrugated 
interface will be derived, and will be solved to first order in the 
corrugation amplitude over the periodicity. Consider first a flat 
interface, where the dielectric constant is for X > 0, Cg X < 0. 

If an incident field is to be totally transmitted as ?E - 0 in the two 
regions, the fields can be expressed in the fonn for p-polarized waves 

A A ikA»/iTsin0,Z ikAV^cose.X 

E = Eq(X sin6j-Z cosGj) e ^ \ > e ^ ^ / X>0 

a ikA»^sin0,Z ikAv^cosO.X 

H E„y e ° ^ ' e ° ^ * (IV. 3.1) 

(i) u 




IkQi^^sInOgZ -ikQ^^cosOgX 


V * E^(X slnOg+Z coseg) e 

OF* POOR QUALITY , ^ ciy -ik^t^cosO^X 

n ^ tQi e e 


X < 0 
(IV. 3. 2) 


Using Snell ‘s law t^sinQj * i^sinOg and the boundary conditions 
(continuity of the fields in the Y-2 plane) one obtains the relations 


hence 




i^Eq * ^^2^0 ’ "‘^osOjEq - cosS^Eq 


C0S6/ 


* "co s O^ * ^9®^^ using Snell *s Taw, one obtains the relation 


t j ? 

l~ sin'^e, 
12 _ ^2 ^ 

l-sin^Oj 


or sin 0, = 


1 e2+C2 


(IV.3. 


This, as expected, gives Brewster angle constraint for e« and c, > 0. 

2 C l 

If Eg < 0 then sin 0^ is positive but larger than 1 so cosQ is itnaginery 

If the Surface is periodically corrugated, given by X = F(2) where 
F(Z) - F(Z-f^) for any integer n, the Bloch-Floquet theorem implies that 
the fields c8n be expressed as summations over absorbed grating vectors, 
i.e., the grating k-vector given by kg * kQg 


E = I E^(X(y^sin 62 +ng) - Z[e2“(>^sin62‘^^9)^3^ 


ikQ(v€2sinej+ng)Z ikQ[c2-(*^2^'''’®l'^^9)^]^X 


X e 


u r ikQ(^^siri02+ng)Z ikQ[ej-(v^sin0,+ng)^]^X 

H = Y— — e e 


X > 0 


<i) ** "n 

n 


' 7^ ^ Ej|^(X(/e^sin02+ng) + Z[e2-(/G^sin02+ng)^]^) 

i kg( vSgS i n02+ng )Z - i kQ[c2- ( v^s i n02+ng ) ^ X 


V p, ikQ(/c^sin02+ng)Z -ikQ[e2-(^sin02+ng)^]^X 
” ' ^ 0) ^ ^n ® ® 

(IV. 3.4) 
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Here ir n * 0 pertains to the surface plasmon mode n * -1 provides 
a radiated component Into dielectric Cj and hence E.j/Eq is the coupling 
coefficient of interest. 

To proceed we apply the boundary conditions and use Snell's Law 

Define {►^sinO^+ng) * v^sinOj^ - law), 

[ej-(y^sin 0 j+ng)^ 3 ^ « v^cosOj^, and [G2*(^sinfl2+ng)^3^ = v^cos02p. 


If 0(2) ^ the boundary conditions become 

(sinOj^sin0(2)-cos0jj^cosO(z)) 


ink_z ikn^^cosO. F(z) 
Ett Pn « ' « ® ‘ 


^ ink.z -ikn/ipcos09/(z) 

s J Ejl^ e ^ e ^ (sin02pSine(z)+cose2^cos0(z)) , 

ink z ik^v^cos , F(z) ink 2 -ikQV^cosO, F( 

Ht'. I En e 9 e 0 1 in ^ J E;, e 9 ^ ^ ^ 

” " (IV. 3 . 5 ) 

Specializing now to a sinusoidal corrugation F(z) * -I- sin k_z , then 

g 9 

ecosk„z 


cos0(z) --7 — p- — p — r, sinO(z) = -7 — p~* 9 ^ 
/He cos k^z /l+e cos 




(IV. 3 . 6) 


Substituting into the boundar-y condition expressions above and using a 
well known expansion one oboaens 

= J J (a) . one has 

m ^ 


. 1 S * 


n,m 


Z 

9 cos 6 j„) (csinOj^cosk z-cosOi„) 

9 


i(n+m)k 2 c/cZkfs 

^ ^n ® ^ ‘^m^“1c COS02,,) (esin02„coskgZ+cosO2^) , 

^ (IV. 3 . 7 ) 
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This can be expressed, using ^ associating 

wave numbers as 

e/c7ko (N-n)k„ 

J VN-n^"n^ ^7e^kQCos0j^ " 

Ei/Eo^Q (^“^)k- 


(IV. 3. 8) 


and 


p i/p "l( 

I Vn.„(^ ^^Oln) = ^ I E;0N-n(^ 


e/e^ko 


which must be valid for any integer N. 

The System above can be solved in the small c limit. To first order 
in e, Eq.Eq are the same as in the flat interface case, and Ej^ E|, E j, 

E^j, the only other nonvanishing components in this limit, can be determined. 
The small argument of the Bessel functions is used, i.e., Jn(cX) -> 1, 

J^j(ex) , all others vanish to first order. 

The equations for N=0 yield the flat interface relations between 

Eq, E^ and hence yields (the terms involving E^^j, E|j 

may be neglected as they are second order in c). The equations for N = ±1 
are 

Ev'^c k k 

N=l, E^: E^{-cos6jj) + Eq(— cosB^q) (yg^i^^cosO^Q " cosB^q) 

' Ei(cos62l) + "'820) ( ^gkocole;,) ^ , 


N= 


“E./c ”k 

=-l, E^: E_j(-cos6j_j) + ^0^ 2kg ° jkgCosBjQ ‘ ^°sB^q) 


(IV. 3. 9) 
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e:i(cos02.i) + E^(-^ + cosejo) 


c/eTkrt c/^kn 

N«l, H^; »^Ej+Eo(-^cos0jq)] = /^[E'+E^(-^^cos02o)] • 

'9 9 

-c^n “=*^0 

N=-l, H^: /cJ[E_j+EQ(~2j^ cos0jq)] = i/e^E'j+EQ( 2k cosOgg)] 


The radiated component of the field is proportional to E_j. The 

/iT 

equations of interest can be expressed, with Egy'^ Eg, 

F o p'] 9 

»E_jCos0j_j + ^107^ ^^"~k 

9 ^ d g 


G k 

*^E_, " 2^0 ^^"k“ ^ ^-1 

• A y y 


(IV. 3. 10) 


These can be combined to eliminate E^j, yielding 


Crij 10 2 


*^*^0 __2 


cos©2_j kg 
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7 02 / “^2 COSOpg 

Using sine^g =/-^ . cosOjg =/-^ , 


etc., one finds 




(IV. 3. 12) 


and 0j__j and ®2+l provide the angles for the emergence of the modes. 
When U2I » |cjj , this reduces to (for Cj = 1) 
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— , and if « l 
i '^0 
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<< 1, one finds 


-1 


E»e 


1 /d x2 


12 


{hr = cc k, 


'0 


The coupling for a small e can be substantial due to the large magnitude 
of . akQ is a simple approximation of the coupling coefficient deduced 
by Seymour et al^^^, who have increased otko bv increasing c]. 


The above results provide detailed coupling parameters. In particular 
at optical frequencies they indicate that 1% coupling can be achieved 

o 

for a 300 A height in the grooves with a spacing of the order of a wavelength 
The basic problem with this particular coupling scheme is the 
optimal design of the corrugation to prevent scattering into higher 
order modCa, A sinusoidal grating Is in particular not optimal. Seymour, 
Koteles, and Stegemen^^^ have experimentally demonstrated a coupling 
efficiency of 8-95;; using overcoated gratings to help confine the radiation 
field. This effort was at 118 y and does demonstrate possibilities for 
visible frequencies. 

The overall power efficiency for the coupling for the particular 
mode can be estimated from the theory of Tien and Uhlrich^^^^ for integrated 
optical coupling. The coupling efficiency is given by 

l-exp[-(3"+akQ)£] 2 

• ^^0^- • (g'Vak,o)k > (IV. 3. 14) 


where 3'' is the inverse of the plasmon decay length. ,and_g, the grating 
length. At visible frequencies for which 3“ is large n can be s 1005i 
for l > (akQ)"^ however the energy is absorbed. One must guarantee that 
i « 3 '*"^ and akQ » 3" for efficient loss free coupling. This is difficult 
to achieve at room temperature; however, not impossible with careful design. 
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An alternative solution to the coupling of radiation to tunnel 
junctions is the utilization of tapers in analogy to microwave coupling^^®^. 
Due to the slow wave nature of the surface plasmon mode and the abrupt 
ending of conventional tunnel junctions, a large mismatch in wave impedance 
exists between the junction gap mode and the surrounding media. This 
implies that little light can be transmitted or received. By gradually 
taping the ends of the device a wave impedence transformer is made whereby 
the gap mode becomes radiative. 

Letting the height of the gap, Xq be sufficiently small such that 
only the surface plasmon forward and backward modes exist in the junction 
(a good approximation, as Xq s 10 A) a coupled mode analysis' ' can be 
performed to obtain the coupling coefficient, k, between these modes. If 
K is high, as in the conventional junction, most of the forward wave will 
be converted to the backward wave (high reflection) and not much light vtfill 
enter or leave the device. If k is low» coupling between the junction and 
surrounding media is high. 

The transverse electric and magnetic fields of the forward and back- 
ward waves, Eg, Hg and Ej^, Hjj respectively, may be represented by. 


Eg == A(z)eg(x,z) 
Hg = A(z)Rg(x,z) 

^ B(z)e^(x,z) 
Hb = B(z)Bjj(x,z) 


(IV. 4.1) 


(IV. 4. 2) 


where A is the forward modal amplitude coefficient, and eai jia ^Te the 
transverse vector modal functions. B, have corresponding definitions 

for the backward wave. Quantities e, h, are normalized using 


in,n * a,b and S Is the cross-sectional area of the junction. 

Denoting 3 as the modal propagation constant, the coupled mode dif 
ferential equations governing A and B are. 
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^ + iBB = ikg^jA 


(IV. 4. 4) 
(IV. 4. 5) 


where k^j^ is the coupling coefficient from A to B, and kj^^^ the same from 
B to A. By assumption, only two modes exist in the junction, therefore 
their combined power is conserved. Consequently, k^jj " *^ba ‘ * 

In tapered waveguides the guide medium changes as a function of 
propagation direction. The varying boundary conditions allow the normally non- 
interacting (different spatial frequency components) forward and backward 
modes to couple. As a result, 3, e, R, all become functions of the pro- 
pagation direction, z. The coupling coefficient is shown to be given by^^®^ 



> V <*= = Kab 
s 


(IV. 4. 6) 


The field solutions for the antisymmetric plasmon mode given in 
Appendix II are good approximations for the vector modal functions. The 
limitations are Imposed by the boundary conditions taking into account 
the variation in Xq(z) with z. Considering continuity of tangential H 
and E, 


[H ].? = 0 
z Xq 


[E 


dXr 


1 


X dz we 6x ■'x 


m.1 


0 


(IV. 4. 7) 
(IV. 4. 8) 


Eqn. 7 is as before. Eqn. 8 contains an additional term due to the sloping 
bbundary x«. If this term is small compared with the second, i.e., if 

dXQ j 

<< ^ ' f" «'hen the boundary conditions are identical to those of 
section (IH-ii;. As a consequence, the aforementioned fields become 
good approximations to the transverse vector modal functions. Using, the 
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fields given In Appendix II , away from resonance where C £ 1, we have 


and 


h “ cosh(v^x-*„)S 

1. ■ 

"b ■ coshCv^x . 



tanh 
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4>, 


m 


tanh 




(IV, 4. 9) 
(IV. 4.10) 


Comparing the power in the insulator to the total power for the 
visible region of the spectrunu it is noted that the fields are concentrated 
In the insulator (C s i). Therefore, essentially all coupling will occur 
there, and only the fields in the insulator are used to calculate the 
coupling coefficient. Using Eqns. (6,9 and 10), assuming small variations 
of the fields across the barrier. 


dVfa (2Vf,Xg-3(|)g)xQ 

= 'ir 6 


(IV. 4. 11) 


For the visible region of the spectrum, Vj^ may be approximated by, 




(A 


+(p ) 


(IV. 4. 12) 


'^*soo» '^s* ‘^’m ^ ^ respectively. 

Using this expression in Eqn. (11), one obtains 


K = 




6x, 


dz 


(IV. 4. 13) 


‘I’ 


S" 


For a symmetrical structure = t|)„_ and v^. = - — . Thus 

y 


- 99 - 


(IV. 4. 14) 




!!p o 

^ * 0 ^ 


For an exponential taper with profile Xq 
value of K of 

2 

l*=pe«kl 


one obtains a peak 


with a In cm'*. Here we have taken Wp/c s 5x10"' cm * and A s 100 A being 

the maximum value of Xn for which the approximations are valid. Thus 

-1 ^ 4-1 

Kpgak « 4a s 400 cm for a = 10 Cm . This indicates that an exponential 
taper of the order of 3 pm in length will guarantee low power coupling 

o o 

into the backward mode while going from 100 A to 10 A. More detailed 
numerical calculations of the power converted into the backward mode are 
presently being carried out. 


and A s 100 A being 
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IV-5 Heterodyne Diagnostic Experiment - Numerical Estimates 

To investigate the non-linearity of tunneling junctions at high 
frequencies as well as the various coupling schemes discussed in the 
previous sections we have devised a heterodyne mixing experiment which 
uses conveniently available sources in our lab, and provides a high fre- 
quency r.f, output which can be monitored on a spectrum analyzer. This 
involves the heterodyne mixing of the output of a water vapor laser with 
frequency (llTHz) and the output of a carbon dioxide laser with a freauency 
of 32 THz. The output frequency wjp ~ should be at 

approximately 20 GHz (Fig. IV-8). The geometry and fArication of the 
devices is discussed in section V-5. The preliminary calculations are 
carried out in three stages, First the input coupling is analyzed, then 
the conversion efficiency, and finally the output coupling to the spectrum 
analyzer. 

The first step in analyzing the input coupling is to determine the 
incident power density, We have assumed 75 mwatts of water vapor laser 
power, and 1.2 watts of COg laser power available. Both of these levels 
have been achieved in our laboratory. We believe that the water vapor laser 
power can be increased, and that 75 mwatts is a conservative estimate. 

The 28 pm(llTHz) beam is to be at almost normal incidence to the silicon 
substrate and will therefore suffer a reflection loss p 


P 


2 





(IV, 5.1) 


of .274. The 9.3 pm (32THz) beam will be incident at an angle of 40" from the 
normal , the difference in angle being determined by the relative 
locations of the gain maximum of the antenna at the two frequencies. 

At 40® the reflection loss is 



COS02 -ng^cos02 
COS02 +n^^.cos02 



.210 


where the angles are shown in Fig. (IV-9). 


(IV. 5. 2) 
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IV-8 Schematic diagram of the heterodyne diagnostics experiment. 
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IV-9 Directions of first two lobes in the 9,3 pm antenna pattern. 
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The diffraction limited spot size can be determined from the beam 
diameter (D), the wavelength, and the focal length of the lens or mirror 
being used for focussing. The radius of the spot Is approximately 

rp = 1,22 ^/f^+(|)^ (IV. 5. 3) 


which Is 114 pm for the COg laser and 683 pm for the KgO laser, The 
available Incident power density at the antenna Is then * ^laser 

|1“P^| X CirrQ]’^ = 3.7 watts/cm^ HgO power and 2.4x10^ watts/cm^ COg power. 

In order to estimate the fraction of Incident laser Intensity coupled 
into: the antenna structure and hence coupled to the junction, we next 
determine V^, the voltage Induced across the antenna, which requires 
estimates of the antenna gain at the frequencies of Interest and for the 
particular antenna lobes of Interest. 

Consideration of the true dielectric constant of the metals yields 

» 

Eqn. (IV-'15) as the anal iKgiye of a low frequency antenna calculation. The 

2 

magnitude of the first order corrections to Is of the order of (w/w_) 

^ p ^ ' 

X k«. For nickel , the metal currently in use, » 3x10’' at 

I .4 P 

28; pm and 2.4x10 at 9,3 pm thus k ~ kp. Also, the decay length 
(the length for which at « 1) Is in the centimeter range for both of these 
wavelengths. Since all of the antenna structures are much smaller than 
this, at can be assumed to be much less than 1. 

With 0/k s 1, and at « 1 for both frequencies of Interest, Eq. 

IV-15 reduces to the normal low frequency antenna result. This Implies 
that standard antenna theory can be used to calculate the gain of the 
antenna structures used in the length confinement coupling scheme. We 
have carried out a computer calculation of the gain at 28 pm [18 
The results indicate that for an antenna of length L “ 5.5A, the gain 
should be approximately 34 at an angle which is nearly normal to the pre- 
ferentially etched edge. At 9.3 pm, formula IV-15 can be Integrated 
over solid angle to give the total power radiated, cos(kt) = 0 for 
antenna of length kQl = 337i so that 


P = J1 
« 8 , 


2ir 7T /cos(kf.LcosO)fl\^ 

/ ^ / [ — hw-' J 


(IV.5.4) 
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substituting r ■ cos6 gives 
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P„^ 



(IV. 5.4) 


Splitting— into - fcppy 2T[TTrT > substituting v« 1 - r 


and u « 1 4- r Into the subsequent Integrals gives 



,2|. ? 2 l-cos33irv ^ . ? l-c os66iiv 

O'- ^ V ^ ^V) 


d\i] 


■ ^ i^[2nn667r+.5772-C1(6670) - K In(132ijf.5772-Ci(132ii)] 

nlfl 

- (IV. 5. 5) 

The gain at any angle, Is then 


(Pp/4tT ” 8.52 


cos 


(KqUosGI+I 2 
^ 'i ' 5 


sinO 


(IV. 5. 6) 


Maxima in the pattern are near angles for which JS.57icosO - niir, 
m = 1,2,... il6. Letting m = 15 yields for the gain of the second lobe 
5.44, This gain calculation is for an antenna in a uniform surrounding 
medium. In the actual structure, the antenna current is at the silicon- 
metal interface, and, according to the results of section (2) radiates 
very little into the air half space. The total power radiated should 
thus be half of the above doubling the gain to 10.88 at an angle of 
24.4°. The radiation strikes the etched Si edge at an angle of 10.7 
degrees (Fig. IV-9) with respect to the normal. Snell’s law, shows that 
the incident angle of light required to couple into this lobe is 40° (op. 

It is not possible to couple to the primary lobe since the critical angle 
Is exceeded. 

At both 9.3 microns and 28 microns, the response of the metal is large 
enough and the junction small enough that the surface plasmon wave gap mode 
reduces to the quasi-static limit of a lumped capacitor. The result is 
that the total diode Impedance is 
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where C = 


> sA 


and Rq Is the dynamic tunneling resistance of the diode. 

Assuming A ^ 1 umx200A, d ® 20A and c^ef 2 yields C = 1,77x10"^® farads 
whdre c^gj^ is an estimate of the high frequency dielectric constant of the 20 A 
nickel oxide. These formulas lead to the following impedances for the 
diode at 28 urn and 9.3 pm. 


Zp (2 « [13.82 - j81.96]n 

Zp 19.3) * [1.59 - j28.15]R 


(IV.5.8) 
(IV. 5.9) 


where Rp = 500fl. 

Integral wavelength antennas have current nodes at the feed point. 

This leads to a large real antenna impedance (500-2000R). Shortening 
the length of the entenna slightly results in a slightly inductive impedance. 
For the worst case calculation, we use a 2000f2 impedance and maxiinize 
5ntenna-diode mismatch. For the best case calculation we assume the 
<sntenna to be slightly inductive to resonate out the capacitive 
parts of the terminations. 

Standard antenna theory gives^^^^ 

2 

(IV. 5. 10) 


V^ 

p - rr - 1 A 
‘5!F®’inc 'Sr 

a 


for the maximum power coupled tp the load. This assumes a load conjugately 
matched to the antenna. From this formula, one can solve for V^ and 

a 

calculate the actual coupled power as 




= ( and 

I'', I I^dI 

■ 


(IV. 5. 11) 
(IV. 5. 12) 


(see Fig. IV-10) 
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IV-10 Equivalent circuit for antenna coupling to the load and the I.F. 
connection. 
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The above numbers yield for the maximum and minimum induced voltages. 


Maximum 

Vp(C02) = 148 mv 
Vjj(H 20) - 90 mv 


Minimum 

Vp(C02) = 52 mv 
Vjj(H 20) = 46 mv 


Knowing the load (diode) voltages allows us to estimate the inter- 
mediate frequency (19,56 GHz) signal generated by the diode. 

Using a Taylor's series expansion about the d.c. bias gives 

(IV.5.13) 
n "■ 

where V is the total junction voltage and 


i(") 


3V 


The n=4 term yields ^ I 


(4) 


(IV. 5. 14) 


a 3 ^ Sw. *u)p ) t 

V contains the term 4V ((jj^)V{( 0 o)e * ^ which contributes the 

20 GHz term. As a minimum 1^^\ we take Javan's^^^^ number of 10"^ tt 

XO M!Ei_ ^ Our devices are noticibly more nonlinear (Fig. IV-11), 
(mv) 

the sensitivity S = ranging from two to six per volt. Those 

of Javan are in the range of 10"^ per volt. Numerical calculation 
gives a minimum of 2x10 amps (mv) for the of our planar diode, 

using a low frequency I-V curve. There is some 

experimental evidence^^^^ which suggests that the planar diodes are more 
nonlinear at higher frequencies. In addition, we hope to be able to 
increase the device nonlinearity, as already discussed. 

Using 3x10”^^ ampy(mv)^ as a reasonable maximum for I^^^ results in 
^max^^®'® GHz) s 1x10’^ amps and I^.^(19.5 GHz) s0.7xl0"^ amps. The 
current is split between the diode, and load (Fig. IV-10) as 

T - T w 500^ 

Moad ■ MF’^ 550n 


I 

Y 
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IV-n Planar junction current-voltage characteristics. 


and so the 20 GHz power launched down the coplanar waveguide is I 

f 

^TL "" ^load ^ ^ 

Thus for the maximun and minimum values of one obtains 



-4 

= 1.15x10 mwatts = -39 dbm max 
*8 

= 2.8x10 mwatts = -75 drub min 


There are several additional sources of loss in the diode to 
spectrum analyser I. F. connection. These include, propagation loss in 
the co-planar waveguide (CPW) taper, propagation loss in the coaxial 
Cable which connects the CPW to the spectrum analyzer, and losses at 
each of the microwave connectors. Using the supplied information on 
the cables and connectors we have estimated these additional losses 
to total less than 7.1 db. Subtracting this from the IF power level 
sent down the transmission guide yields the final minimum and maxi- 
mum expected IF signals at the spectrum analyzer as Pg^(19.5 GHz) = 

-46 dbm maximum and -83 dbm minimum. 

The Tektronix spectrum analyzer which we have will detect signals 
as low as -95 dbm , if the bandwidth is sufficiently narrow. At its 
widest bandwidth the sensitivity is -65 dbm. Phase sensitive post detection 
can also be used to increase the sensitivity further. 

These estimates certainly indicate that this particular experiment 
should be succesful with out continuous wave water vapor and CO 2 laser 
sources and our present devices shown in Fiq. (IV-12). 

Since the actual structure consists of an array of two diodes 
raised from the substrate, the coupling characteristics will differ from 
those calculated above. The estimates are not expected to change sig- 
nificantly due to the distance between the substrate and the center electrode 

O 

since this is of the order of a few hundred A, which is appreciably less 
than the plasmon wavelength. The thickness of the center electrode is 
however less than the penetration depth of the single surface plasmon 
or of the same order and thus the mode, and hence characteristic impedance, 
will be altered somewhat. In addition, the short stub (center electrode) 
causes a phase shift so that the junctions are shifted from an end 
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IV-12 Photographs c/ metal -barrier-metal tunneling devices fabricated 
at Berkeley, a) Overall view, b) Magnified to show antenna, 
c) Further magnified to show the junction. Antenna width is 
1 pm. -Ill- 




termination. As long as the resonance condition is satisfied, this 
should not appreciably alter the expected signal strength. We are 
presently analyzing these two effects. 

Finally, shifting the experiment towards optical frequencies is 
expected to change the results in several ways. The devices and antenna 
structures will have to be made smaller or traveling wave junctions will 
have to be exploited. Since the plasmon becomes more lossy as one enters 
the relaxation region or in the visible is even influenced by the resonance 
behaviour of the plasmon, the coupling will decrease. A great deal could 
be learned by considering the coupling as a function of temperature at 
these short wavelengths. However, it is expected that no fundamental 
problems exist for coupling at visible frequencies and we anticipate 
moving in this direction when our basic results have experimental veri- 
fication. 
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V FABRICATION OF JUNCTION DEVICES 


V-1 Introduction 

In the present chapter we discuss the fabrication of junction devices 
presently being used to investigate various coupling and tunneling aspects 
pertaining to direct conversion through tunneling. 

We have generally been using three configurations for our tunneling 
devices: the metal-oxide-p semiconductor, the metal -oxide-ptype semi- 
conductor-oxide-metal junction device, and the metal -barrier-metal 
junction device. The first two have been useful because of their inherently 
fast response time and thus for the detection of short optical pulses. 

This requires no coupling to the surface mode and is therefore much easier 
to demonstrate than a highly efficient direct conversion process. Any 
utilization of these junctions to exhibit conversion efficiencies theore- 
tically available from tunneling devices will require the development of 
such coupling schemes. 

The primary focus of the present chapter will be the edge metal-oxide- 
metal device presently being utilized In a fourth order mixing experiment 
being implemented as a diagnostic tool to measure both the junction non- 
linearity and coupling at high frequencies. Unfortunately the metal- 
oxide-metal (MOM) is not a highly nonlinear junction, but it is expected 
that it can be improved as fabrication procedures evolve and possibly 
by combining it with semiconductor technology. Prior to describing its 
fabrication in detail, we present a short summary of the work on the 
two other devices. 

V- 2 The M-O-P-O-M Junction 

The M-O-P-O-M (Fig. (V-1)) interdigitated electrode structure presently 
operates well as a fast conventional photodetector with a response time 
of approximately 50 picoseconds^^^. The carriers are generated in the bulk 
semiconductor by flooding the junction with the incident radiation through 
the interdigitation. This interdigitation can in principle be used as a 
grating coupling structure for the excitation of surface waves however 
this is not essential for the detector, hence its simplicity, although we 
have still referred to it as a surface detector. This is the simplest of 
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-0-M device with interdigital electrode 




the three devices to fabricate and will conclevably be utilized as a 
surface device for possible conversion experiments by attempting this 
coupling into the oxide layer. 

To fabricate the device the silicon wafer is first cleaned and a thin 
oxide layer is grown thermally in a dry oxygen atmosphere. The thin oxide 
thickness is monitored using an ellipsometer, typical thicknesses presently 
ranging from thirty to sixty A. 

Subsequent to oxidation » both aluminum electrodes are evaporated and 
patterned using a lift-off technique. This particular tunneling device 
operates by having electrons tunnel from one eleptrode through the oxide 
into the semiconductor substrate. In order to complete the electron path 
from one electrode through the oxides, to the second electrode, free 
carriers must be supplied to the semiconductor either by excitation or by 
doping. For the detector, we utilize photo-generation. For a tunneling 
converter this would be a heavily doped layer or even conclevably a metallic 

layer. We are attempting to modify this structure by first laying the 

0 

device on a saphire substrate with approximately 3000 A of silicon. Such 
a thin Si layer would confine the current path. Requiring no new fabrica- 
tion steps, the simplicity of the device is retained, while allowing the 
eventual fabrication of these devices on a metallic substrate for the guiding 
and coupling. 

V-3 The M-O-S (p'*'^) Junction 

There are two important factors motivating the fabrication of the 
second class of (M-0-S(p'^'*^)) devices. These are expected to have more 
nonlinear current voltage characteristics than simple M-O-M devices and 
hence provide better conversion characteristics. Secondly the plasma 
resonance frequency can be tailm’ed by doping, thus possibly providing a 
means of enhancing conversion characteristics. Three types of semiconductors 
(Si, 6aAs, SnTe) are being used. For silicon, thermal oxidation is first 
carried out. Subsequently the oxide is patterned, and the silicon ion 
implanted to create a p^^ region. The wafers are then annealed, and a 
thin oxide layer grown and measured as for the m-o-P-O-M devices. The 
metal can then either be evaporated or sputtered to form the electrode 
structure. If p^''^ 6aAs is used, it is obtained either by purchasing it 
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directly from a supplier or by growing epitaxially. SiOg pr SigN^ is then 
deposited and a window opened through the insulator (Fig. (V-2)). Ga^O^ 
is then grown using the sputter oxidation technique described below for 
the "edge MOM" fabrication procudure. The metal layer un then also be 
deposited in the sputtering machine without breaking the vacuum. 

SnTe has the convenient property that it can be evaporated and still 
retain semiconducting characteristics, the resultant film being 
This allows the deposition and patterning of the metal layer first and 
the use of a metal oxide tunneling layer obtained via sputter oxidation, 
thermal oxidation in air, or by placing it in an oxygen plasma in an 
evaporator. The last method allows the evaporation of SnTe onto the 
substrate, or immediately following the oxidation. To finish the device, 
the SnTe must be patterned, either by lift-off or chemically using HCf . 

V-4 The M-O-M Tunneling Junction Fabrication 

The EDGE MOM is so named because the metal -oxide -metal junction is 
forriied on the edge of one of the two metallic electrodes. In this con- 
figuration, one dimension of the device is thus determined by the thickness 

of the metal on whose edge the junction is fabricated (Fig. (V-3)). This 

0 . ' 

dimension may be as small as 100-200 A. This device is fabricated using 
two rf sputtering deposition steps. The first electrode is patterned 
using an ion milling machine, and the second with a wet chemical etchant. 

In our fabrication procedure for these junctions rf sputtering 
has been heavily used. Three different processes are involved; sputter 
etching, gputter deposition, and sputter-bias deposition. All three of 
these modes involve a plasma whose ions are accelerated into either the 
material to be deposited (target), or the substrate, or both. In sputter 
deposition, the ions are accelerated into the target where they dislodge 
target molecules which drift down onto the substrate. 

When sputtering a metal, a purely nonreactive Ar plasma is used. When 
depositing SiOg (or Si^N^), some Og (or Ng) must be added to the plasma 
tp keep silicon compounds with fewer oxygen (nitrogen) atoms from forming. 
Iii sputter etching, the plasma ions are accelerated into the substrate, 
with resultant ballistic etching. Ions are accelerated into both the 
target and substrate in the sputter bias mode. This deposits and etches 


I 

f 




- 118 - 


ORlG^NAt 
OF POOR QUAUTY 


\r Metd 


Si,N^ 

V/ 

Si,N4 


^Go Oxide 

GoAs - 



M-O-P^^ (GaAs) structure used to study interactions between surface 
electromagnetic waves and tunneling current. 
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V-3 Edge-MOM structure, a) Front view b) AnguUr view 



simultaneously resulting in smoother, lower resistance metal films. 

The sputter deposition mode is used to clean each target before the 
fabrication process. Material is sputtered off the target onto a blank. 

Any absorbed impurities are then sputtered away before rotating the sub- 
strate under the target. 

The first electrode consists of an insulator-metal -insulator sandwich, i 
The substrate is introduced into the sputtering machine, the targets are 
cleaned, and the insulator is sputter deposited. If SiOg is used as the 
insulator, the metal target to be used for the electrode material must 
be recleaned to remove any nickel oxide which grew in the argon-oxygen plasma 
A thin layer of metal is deposited using the sputter bias mode. It is 
important to obtain the best possible metal film during this step, since 

O 

it is so thin, typically 200-400 A. ^he upper insulating layer is then 
deposited. The lower insulating film protects the srTIcon surface during 
a subsequent ion milling procedure, and the milling process is easily stopped 
before etching all the way through this layer and into the wafer. The upper 
insulating layer keeps the two electrodes from shorting. This insulator 
must be made thick enough to minimize the electrode to electrode capacitance, 

o 

typically greater than 1000 A. Photo resist is next spun on and patterned 
by standard photolithographic procedures. The first electrode is then 
ready to be etched. 

The sandwich structure is etched in an ion mill, in order to obtain 
smooth vertical walls. This is essential to building devices on the edge 
of tfie metal. The photoresist used to pattern this electrode must be post- 
baked at 120° for at least 15 minutes, to harden it prior to being placed 
in the ion mill. The Si wafer is also heat sunk with gallium to a 
cooled substrate table. If the substrate is not heat sunk, two problems 
arise in the etching. First, the photoresist becomes too hot, and 
polymerizes, making it very difficult to remove. Secondly, the etching 
process is disturbed. The part of the wafer being etched, develops "hot'' 
and "cold" spots which tend to either etch quickly, or very slowly. This 
results in parts of the wafer effectively not being etched, while others 
are etched too deeply. The etch process is timed to etch into the lower 
insulator layer, but not through it. Etching through the lower insulating 
layer results in a roughening of the surface of the Si, making subsequent 
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chemical etching of the upper electrode difficult. This problem cannot 
be eliminated by lowering the Incident Ion energy, as even low energy ion 
etching has been observed to cause surface damage^^^ 

At the end of a usual etching process, the photoresist can be removed 
by ultrasonically cleaning the wafer in acetone. Occasionally, it Is also 
necessary to do a short plasma ashing to remove the last photoresist 
residue. The wafer is now ready for the upper electrode deposition. 

Subsequent to the lower electrode patterning, the wafer is returned 
to the rf sputtering system. Without breaking the vacuum , the wafer is 
^putter cleaned, the oxide barrier grown on the lower electrode edge, and 
the upper electrode deposited. The wafer is sputter cleaned (etched) for 
two minutes in an argon atmosphere to remove any oxide which may already 
have grown on the lower electrode. This also removes a small amount of 
residual photo resist if the stripping process has not entirely removed 
it. The lower electrode is then sputter oxidized, and the upper 
electrode metal deposited. 

The rf sputter oxidation procedure developed by Greiner^ involves 
two competing processes. The substrate is sputter etched using a plasma 
containing oxygen. The etch process removes material from the substrate at a 


constant rate R. If X is the thickness of oxide present, = -R for the etch 
process alone. The presence of oxygen in the plasma causes metal oxide 
to grow. However, the rate of growth is not constant, since it requires 
metallic atoms to be available at the growth surface. This exposure 
decreases with increasing oxide thickness, and the growth rate decreases. 

The oxide growtjj^jjate in a plasma decreases as — = Ke'^^^0, 

giving = Ke ^ - R for the net rate of change in oxide thickness 
during sputter oxidation. As metal oxide grows, the growth rate slows, 
until it equals the etch rate, and no further oxidation occurs. Setting 
3X/3t = 0 gives = XqIo^ as the equilibrium 

thickness. This thickness is approached exponentially with a time 

: ' . i 

constant of Xq/R. This time constant varies substantially from metal 
to metal. In some metals, the time constant, or equilibrium thickness 
are such that for the fabrication of practical tunneling barriers the 
equilibrium is never reached and the time of growth becomes an 


Important parameter determining the barrier thickness. In the case of 
Ni, the constant has been determined to be only a few mlnutes^®^, so that 
the equilibrium thickness is reproducibly achieved by allowing the process 
to run for 10-15 minutes. The equilibrium thickness can be varied by adjusting 
the rf power applied to the substra_te, and the amount of oxygen in the 
plasma. With Nii the equilibrium thickness can be adjusted to produce 
edge MOM's with zero bias dynamic resistances from about 750 q to about 
one Mil. Wafers with lower zero bias resistances have a noticably larger per- 
centage of shorted devices on the wafer. 

The oxide chamber is covered as quickly as possible after it is grown. 

The chamber is evacuated to remove the Ogt and refilled with a pure argon 
sputtering atmosphere. Ni is first sputtered off of the target, onto a 
blank table, to remove any nickel oxide formed on the target. A thin layer 

O 

of Ni (approxi,.'’ately 200-400 A) is then sputter deposited very quickly onto the 
substrate ( approximately two minutes). After this initial layer, which 
protects the barrier, the rest of the upper metal layer is deposited using 
the sputter-bias mode. 

The upper electrode is etched using a standard chemical etchant 
rather than ion milling. It in fact has been experimentally determined 
that ion milling of the second electrode causes the junctions to short, 
it has also been observed that if during the ion milling of the lower 
electrode, the lower insulating layer is etched through, the Ni deposited 
for the second electrode adheres too strongly to the roughened Si sur- 
face. The roughened surface makes it difficult for the chemical etchant 
to thoroughly remove the Ni. To do so would require the wafer to be 
left in the etchant so long that all lines less than 5-10 pm's wide are 
undercut, and etched away. One solution to this problem has been to put 
the wafer in the ion milling machine for a very short period, to etch 
this last layer of Ni. 

With the completion of this final etching, the edge MOM is ready to 
be tested on the curve tracer. This basic procedure for the M-O-M fabri- 
cation has had several and various modifications, and variations which 
have been, or are to be, tried to increase reliability. 

Because of the inherent geometry chosen, (edge), it is easier to 
actually fabricate a two junction series array, than to fabricate a single 
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junction. The steps necessary to contact the buried electrode can be 
eliminated if two junctions are made in series (Fig. V-4), The series 
array, also forms a natural pillbox structure, which can be designed to 
be resonant at the frequencies of interest. For the 28 pm mixing experiment 
two element arrays (Fig. V-4) have been built, with spacings of L = 2BiJn/n 

between elements. These boxes should be resonant, at 9.3 pm as well 
as 28, within the precision of the lithography alignment. Series 
arrays should improve the overall conversion efficiency of the M0M‘s. 

V-5 Fabrication Procedure for Measuring of the High Frequency Nonlinearity 

The procedure described above is for making a single edge MOM diode. 

The current experiment for measuring the high frequency nonlinearity imposes 
additional constraints. The 9.3 pm mixing experiment requires efficient 
input coupling of two laser frequencies, as well as efficient output coupling 
out of the 20 GHz IF. It is also desirable to be able to construct a 
variety of combinations of diodes, including arrays, and balanced mixers. 

The coupling problems and the desireability of multiple junctions has led 
to the incorporation of a variety of additional processing steps, and to 
a search to modify somewhat the basic procedure. 

The coupling considerations have had three main effects on the con- 
struction of the edge MOM devices. The electrode patterns must be 
designed to maximize the coupling. The devices are built near a preferen- 
tially etched edge of the Si wafer and the SiOg is being replaced as the 
insulating layer. These three steps should allow for optimum IF signal 
generation. 

The shaping of the electrodes for good input coupling has been the 
most difficult design problem. Standard antenna theory can be used to 
predict receiving patterns at longer wavelengths. At 28 pm and 9.3 prn 
two problems with this approach begin to manifest themselves. At these 
shorter wavelengths, the metal does not exhibit ideal loss-less free elec- 
tron behaviour. Both the displacement current, and the conduction current 
must be taken into account in calculating any radiation, or receiving 
patterns, as has already been discussed. The limits of optical lithography 
are such that it is not possible to make '-thin wire" antennas without 
using e-beam lithography. While theoretical efforts are being made to 
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V-4 Edge-MOM junction structure showing the buried electrode. 
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understand these problems » there are, as of this writing, no sure directions 
for proper design. 

This lack of a clear theoretical direction has led us to continue to 
build long wire antennas (Fig. V-5). At longer wavelengths, where analysis 
is easier, these show very high gain. There is a need to experiment with 
other shapes, suggested in another section. These antenna structures 
must be incorporated inside a 20 GHz IF electrode pattern, which does 
not radiate. 

The larger electrode structure is that of a coplanar microwave waveouide 
(Fig. V-6) This guide has been designed to have an impedance of 

50 f2's. This waveguide structure offers the possibility of incorporating 
matching elements directly into the waveguide, to maximize the 20 GHz 
coupling. This coplanar waveguide will connect to a taper section built 
on an alumina substrate. The taper changes the dimensions of the waveguide, 
but not its impedance (Fig. V-7). 

This coplanar waveguide structure should couple the 20 GHz IF quite 
readily into the spectrum analyzer. Input coupling considerations, 
combined with the desi reability of a completely integrated structure have 
led to the building of these devices above the preferentially etched edge 
of a Si wafer (Fig. V-8). Calculations show that an integrated antenna 
receives radiation very poorly from the half plane above the dielectric 
(air side). In past work a superstrata, Si prism, was placed on top of 
the antenna. The antenna above a preferential ly etched edge allows the 
substrate to perform the role of the prism and allows one to irradiate 
the antenna at angles which would be totally neglected from the back of 
the substrate. The antenna length is chosen to place the beam maximum 
perpendicular to the preferentially etched edge. The edge is etched 
prior to the device construction. 

The preferential etching of Si is performed by a standard EDP 
(Ethelyne Diamine Pyrocatachol ) process The cleaned Si wafer is 
oxidized in a wet O 2 atmosphere. The atmosphere is changed to N 2 after about 
8000 A of oxide has grown. The oxide is annealed in the Ng atmosphere 
for a half an hour, then removed from the oven. The wafer is patterned 
with photoresist. The SiOg is chemically etched inHF, creating the mask 
for the Si etch step. The photoresist is removed, and the wafer is 
immersed in a heated solution of EDP. The wafer is left in the EDP until 
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V-5 Top view of series array of two tunneling junctions showing two 
4.5 X long wire antennas on either side of the pilT box structure 
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V-6 Three possible configurations for interconnecting the tunneling 
junctions and the microwave network for the heterodyne diagnostic 
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V-7 Overall view of coupling scheme between the microwave spectrum 
analyzer and the tunneling diodes. 
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Diodes above a preferentially etched edge of a silicon wafer. 
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It has been etched through (about four hours). The remaining SiO^ is 
now removed with a buffered HF solution. The substrate is then ready for 
the ROM device to be built on top of it. Having examined the motivation 
for the electrode shape, and the preferential etching of the substrate, there 
remains the choice of substrate and insulator best suited for the 28 pm 


and the 9.3 pm coupling. 

There are only two substrate materials which have been found to have 
the appropriate hardness and solubility requirements for IC processing, 
while also being transparent to 28 pm and 9.3 pm radiation. These are 
high resistivity Si, and high resistivity Carbon (diamond). Diamond is 
both expensive and difficult to preferentially etch or lap. Silicon is, 
therefore, the substrate of choice to Investigate the junction devices. 

We are using > 700 fi-cm n-type Si. It is also necessary to choose an 
insulator which is transparent at 28pm and 9.3 pm. SiOg, the insulator 
we have usually used, absorbs 28 pm light too strongly to be a good choice. 
Experiments are currently underway to try to eliminate SiOg as the 
Insulator. 

The most important place to eliminate SiOg, would be under the antenna 
structure, since this has a strong adverse effect on both the antenna 
pattern and the effective area, Eliminating the SiOg layer entirely is 
not feasible. Although the silicon has sufficiently high resistivity 
tO: isolate the electrodes, one runs into the previously mentioned 
surface roughness problem resulting from the ion milling step. A possible 
way around this is to put down the SiOg layer, and then remove it after 
ion milling. The removal can either be done via chemical (HF) or plasma 


etching. Plasma etching should be done in a barrel reactor, so 
that the ions are nc^ accelerated into the substrate. These procedures 
can be eliminated completely if a more suitable insulator is chosen. A 


suitable material is SuN., which is 
COg and HgO laser wavelengths' . 


reasonably transparent at both the 
The sputtering of SijN^ in an r.f. 


sputtering machine is currently being parametrize^. 


The resultant junction devices are then designed for a high frequency 
measurement of the nonlinear mixing characteristics with an I.F. frequency 
of 20 GHz and an upper frequency corresponding to 10 pm. Various modifi- 
cations to increase nonlinearity and coupling are expected as a result of 
these measurements. 


I 
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' VI OUTLOOK - ABILITY TO REACH THE THEORETICAL LIMIT 


In the present report we have detailed our efforts Involved In an 
investigation of the possibility of utilizing tunneling junction for the 
direct conversion of the free space radiation field to electron flow, 

This has involved the theoretical aspects of tunneling using a stationary 
stpte theory. This has allowed us to provide a unification of the Tien- 
Gofdon model with other models dealing directly with the Interaction 
within the tunneling region. For practical convertors it has made clear 
the analogous pictures of multi-photon absorption ahd Fermi level modulation 
and has provided us with a quantum mechanical view of the Ideal 40^. conversion 
from low frequency rectification. We have seen that the competing 
processes of multi-photon absorption from the field and dielectric absorp- 
tion within the bounding materials Is a problem which must be dealt with. 
Experiments on cooled junctions should help clarify some of these aspects. 

One of the basic problems with regards to the utilization of tunneling 
junctions is the fabrication of a diode which does not have appreciable 
backflow; that is when reverse biased It does not conduct. This seems 
to run counter to the tunneling phenomena itself and the only way this 
appears definitively acheivable at the moment Is by using the energy 
gap of a semiconductor as postulated. The current voltage relationship 
for such a junction should be calculated theoretically to determine whether 
a sharp enough turn-on of the current can be acheived while retaining a 
low back flow. There are also some indications that metal -barrier-metal 
junctions themselves can exhibit a sharp enough turn-on to overcome this.’*’ 

The surface modes (and the coupling to them) which propagate in metal- 
barrier-metal structures have been considered In chapters III and IV 
respectively. We have shown that the coupling between an open surface 
mode and the junction mode Is primarily one of fabrication; that theoretically 
close to 100% coupling can be expected for ideal tapers* Coupling 
between the structure modes and the free space radiation field can be 
accomplished through appropriate antenna design as shown. We have also 
mentioned an end coupling scheme which becomes possible at high frequencies 
and which has demonstrated experlmentaTly. Calculations on the 
Tatter should be carried out as an extension to the present work. 

^See C. Slayman, Ph.D. Thesis (U.C. Berkeley, unpublished (1980)), page 30. 
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The principle problem with surface modes is their loss and 
particularly short propagation length at visible wavelengths. This 
was discussed in detail in section (II1-3) and has been verified by 
numerical calculations (II 1-4). The short propagation length, particularly 
for the junction gap mode can be alleviated by cooling if the relaxation 
time is phonon collision dominated. This is true of the metal, which is 
the most important because of its extremely short relaxation time. 

Picosecond spectroscopy of thin metallic films at low temperatures 
would be interesting since it should provide direct evidence that the 
relaxation time could be extended to at least the picosecond range by 
cooling. Coupling to surface plasma waves and either propagating the 
mode or scattering into the free space radiation field could also be 
used to obtain useful data on losses. 

Heavily doped semiconductors should also be investigated since they 
could possibly be exploited to solve the backflow problem. 

In chapter V, we described the processes which we have used to 
fabricate tunneling devices to be used In our experiments. It is likely 
that future development of molecular beam epitaxy, chemical vapor deposition, 
and electron beam deposition will result in highly improved fabrication 
techniques and will allow much more flexibility in the development of 
semiconductor barriers. 

In concluding this report, one could state that the potential for 
highly efficient convertors is theoretically present. It is true that in 
contrast to depletion layer devices which inherently depend upon the 
internal field to separate the carriers, the upper limit to the efficiency 
is one hundred percent. The reason that thus far the converMon 
efficiencies are so low are due primarily to the following: 

a) The coupling to the gap mode from a single surface mode has not 
been definitively demonstrated let alone designed for efficiency, 
although coupling as well as end coupling and coupling through 
roughness have been demonstrated. 

b) The junction current- voltage characteristic has displayed a 
highly symmetrical behaviour about the zero-bias point for 
virtually all tunneling junctions fabricated thus far. Thus the 
back flow has generally been a problem. 


c) The absorption length due to multiphoton absorption through the 
gap nwde is comparable to or greater than the propagation length 
of the mode in the absence of the multi-photon absorption. 
Theoretically,, perhaps the most crucial problem at this point is a 
calculation of the current-voltage characteristics of hypothetical 
junction configurations which could demonstrate more ideal rectification 
behaviour than the metaVbarrier-metal junction. Eqs.{II.5.6 and 11.5.7) 
could then be used to estimate the absorption length and possible 
efficiencies which might be realized. We have made a start by considering 
the metal -barrier-heavily doped semiconductor device; however this needs 
to be extended to include a further metallic boundary on the semiconductor 
side. 


APPENDIX I PROGRAM FOR CALCULATING THE DISPERSION CHARACTERISTICS OF 
THE SURFACE ELECTROMAGNETIC WAVE GAP MODE 


Input parameters are: 

xss - Normalized freqency w/ti-v where o) is the plasma frequency of 

♦1 w 

the heavily doped semiconductor 
e(l), f(l) - dielectric constant of the film is CQ(e(l)+jf(l)) 
es - dielectric constant of film as ® 


xs - 

ws - a)j/(Og = 1 

^m - dielectric constant of film at w ® 


m - 
m ■ 

s - barrier thickness 


The most important output is dm, dp, rk, rkk, xss and s. 
dm - dm = 0 when a solution is obtained, 
dp - n(2Tr) for a solution (phase condition) 
rk - Re e/oi/pQCQ 
rkk - Im 3 /w/mq£q 


Program Disp.f. Notatio n 


'‘o “ ° P® 

= J8878.3 Xq = 3.3383 (jm Wg = 5.6636 x lo’'' 


8 = rk + jrkk 

2 

xss = (u>/li)n) x2 = (xss) 
s(l) = = (e(l)+jf(l)) 

s(2) = ejj = (e(2)+jf(2))ea 
e( 3) = = (e(3)+Jft3))em 

xs = l/Cs«o 


sk(l), sk(2), sk(3) 
change 


' “ ps ''“0 

• V“0 
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X . k e 

** '' n' ■ ■ “ ' ■ 

W^MGq / (D MCq 0 


Controls 

js, < 0 write 

X. yy, e(i), f(i) 

+ 2+2 
Rey I my 




/ V"i 


C #(i) V I 

^ V 


sk(i) < 0 change b(i) -► 7T-b(i) (rk+jrkk)xss = 




k 

p6 


- 0 test for X < 0 

> 0 change b(i) tt - b(i) (rk+jrkk) =(^)(t^) rk+jrkk 

normalize^ ^o k(w) 


'po' 


a(i) = [Amp(k^-w^PE)]/(po) 1 -► barrier 2) s.c. 3) metal 

b(i) = Arc(k^-oi)^yE) 

g(i) = Amp e{i) h(i) = Arc E{i) 


pp' = < e(i) - < E(i) - {< y(i)-< Y(i)) 

bb = (amp e{i)/amp e(1)) x Amp(y(l)/amp y(i)) 

z(J) = 1 (bb)^ + s(bb) X cos(pp) = squared ampl. of ^|y 


y(j) = squared ampl. of 1- 


iiiu 


^1) Y(i 


rq(d') = <(1+- 


Ell 


eTI 


r(j) = <(1 


Ell 




d = 2 a(l) X S Normalized to kQ (dkg) = S 

cb = cos(b(l)) sb - sin(b(l)) dm ampl. cond. dp -► phase 


PRINTOUT 

Z S(l). 1 S(2). y(l). y(2), z(l). z(2) 

+ + + + 
Amp squared sum 
of dnff 


akl, ak2, ak3, bkl» bk2. bk3 
Ini.(Yi.Y2.Y3) Re(Yi*Y2»Y3) 
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prugram Dibp 
dimension sii(3)«zs(3) 

dimension e(5>»f(5)»a<5)» b<S)* g(5)» h<5)i r (51^ rq<5)» 2 < 5 ) / y(t») 

i - 4 -- xss> x 2 f dp» es* em» xsi rks» ws; xmi rkm» wm* 5 

e(n4c(2).e<3)<fU)4f<2)«f(3).rk»rkk,sg 
ui Vt X I y y ) sa * a ( 1 ) « a ( 2 > > a < 3 ) « b < 1 ) * b ( 2 ) » b < 3 ) 
g<l )4 g( 2>4 g(3)« h ( 1) « h (2) » h (3) 4 p« bb» hh 
z ( 1 > 4 z ( 2 ) > z < 3 ) < ij ( 1 > * y ( 2 ) » ij < 3 ) < d I cb « sb « dm« p p 
cb2^ sb2* cb3« sb3» ak i» ak2» ak3« bk 1< bk2« b k3 
rq<l)«rq(2)/r<l)4r(2)ipi<&s^C4po>ei«e2 


double 
double 
double 
doub 1 e 
double 
doub 1 e 
doub 1 e 


prec i sion 
prec ision 
prec ision 
prec ision 
precision 
precision 
prec isxon 
pO*lB078. 3 
pi>=3. 1415926536 
rk«l. 7 
rkk=. 22 f; 

X6B«, 1 
x 2 ~, 01 
ei>^3. 2 
f 1=0. O 
es=l 1 , 1 

xs=l /sqr t < 32. 5) 
ws=l. 

xm=sqr t < 22 . B/32, 5/1. 54) 
U)m=sqrt <25lB. «22. 8/32. 5/1 
em=l . 


rd 
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54) 


s«9. *5. 6635/3. / 10000. 

1 continue 

write <*» 'Center xss» e< 1 >» f ( 1 >/ es» xs< u»s, em» xm» wm< s*'* ♦) ' ) 
read(-e>«^) xss* el> f 1> es* X5< ws/ em< xm« wm> s 
write(«« *) xss / el> fi< es« xs*wsi em> xm> wm» s 
3 continue 

write (*C Center sk ( 1 )» sk <2), sk <3> •'.♦>' ) 
read<-e««) sk ( 1 >< sk (2)> sk (3) 
write(«i«-)6k(l)isk(2)<Bk(3) 

2 continue 

write (■»> 'Center rk and rkk«tr<ts>ib ",$)') 

read(*»*) r k 1 r k k« tr» ts« i s 

rk=rk*xs 5 

r k k=r k k*x ss 

x 2 =xss*xss 

rks=xs*xs 

r km=xmexm 

e(l)=el 

f (|1 )=f 1 

e<' 2 ) =es 1 . -( W 5 *W 5 ) / < x 2 +r k s ) ) 

f (2>=es*( x5#ws*ws)/( x 2 +rk! 3 )/xss 

e(3)=em*( 1. -(wm*wm)/< x 2 +rkm> > 

f (3>=em*< xm*wm*wm)/< x2-*-rkH))/xss 

u»r k*r k-rk k#r k k 

v=* 2 . *rk*rkk 

dp 10 i = 1.3 

x=u--e(i)ex 2 

yysy-f (i)«x 2 

sa=dsqrt < x*x+yy*yy ) 

a ( i > =d sqr t ( sa ) 

b < i )=das in < yy /sa ) 

if < is)78, 79» 79 

78 write<*»*) x 4 y y 1 e < i > / f ( i ) 

79 continue 

iif (sk < i ) ) 7, 9, 8 
9 if<>-) 5, 6 , 6 
5 if (b<i ) )7. 8 » 8 
7 b ( i )=-p i-b ( i ) 
go to 6 



g < i )=d5qrt (e ( i ) «-e( i ( i XH*< i ) ) 

h(i)=:daBin(f(i)/g(i)) 

if <e( i ) )30« 31» 31 

30 h<i)=pi-h(i) 

31 Continue 
10 continue 

do 20 jell 2 

pp«h(k)~h(l)-<b<l()-b<l)) 
csedsin<pp > 
c«dcoB (pp ) 

bb»(g<k)/g(l>)«a<l)/a(k) 
zs ( j )Bbb 

z ( j )«1. -♦•bb#bb4-2. *bb«c 
g < j )«1. -fbb«bb-2, »bb«c 
hh«bb«-ss 
sz=dsqrt < z ( j ) ) 
rq(j>=dasin<hh/sz) 
if ( 1. +bb#c )49, 40» 48 
49 if (hh>46, 47, 47 
4/ rq( j )'=pi-rq< j > 
go to 46 

46 rq< J )=“pi-rq< J ) 

46 continue 

sys=dsqrt < y < j ) ) 
r ( J )*da5in<-hh/sy ) 
if(l.“bb*c> 53,52,52 

53 if <-hh )54, 51, 51 

54 r < J >=-p i~r < j ) 
go to 52 

51 r ( J )=p i-r ( J ) 

52 continue 
20 continue 

d=2, *a< 1 )*s 
cb-c1coB ( b ( 1 ) ) 
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Bb«d sin ( b ( 1 ) ) 

dm=z < 1 >*z (2)*dexp <2. *d*cb )-y < 1 )*y (2) 
dp*rq( 1 )+rq(2)-r < 1 )-r <2>-»d*sb 
rk=rk/xss 

rkk=rkk/xs5 

dp>=dp/pi 

write <*,*) dm, dp, r k, r k k, XS5, s 
if (is)57, 58, 58 

57 write (*,*) z s ( 1 ) , z s (2) , y ( 1 > , y (2) , z < 1 ) , z (2 ) , d 
liirite (#,*) a< 1 ), b ( 1 >, a<2), b(2), a<3), b(3) 
cb2»dcofs(b(2) ) 

Bb2=dsin(b(2) ) 
cb3=dcos(b(3) ) 

»b3=dsin(b <3) ) 

ak l=a ( 1 )#sb 

ak2»a(2)«Bb2 

ak3=a(3)«-sb3 

write(*,*) akl,ak2,ak3 

bkl=a<l)*cb 

bk2=a<2)*cb2 

bk3=a(3)*cb3 

c2=e< 1 )*(e(2)/g<2)/g(2)+rkk*f <2)/g < 2 )/g ( 2 )/rk >/bk 2 /s 

c3=e< 1 >*(e(3)/g(3)/g<3)+rkk*f <3)/g (3)/g (3)/rk >/bk3/s 

c = l. /(I. +c2+c3) 

r kB==r k / (r k*r k+r k k*r k k ) 

rke=rke*e ( 1 ) 

write (»,*) bkl,bk2, bk3 

write<#.*> c,c2,c3,rke 



APPENDIX II THE ANTI -SYMMETRIC GAP MODE 


OR POOR QUAinrY, 


(u) r 

Using Fig. (III-7) and e ' e^( 1 - j £„,,e can represent 

* 10 e 

'•■fn.s 

the fields of the gap mode as 


S(x s |) 


vs(x-f) 


y‘ V 

3Hg vglx-l) 

V5(x-|) 

Z* (OEg ® 

Vg = (a^-Acg)^ 


b(-|sx£^) 


m(xs-^) 


H|^cosh(+v. (x 4)-4 ) H^e ^ ^ 

d 

— cosh(vj,(x+|H5 

D fn j 

— sinh(Vj,(x*f)-*s) -^e 

D m 


V(, = 




Continuity of Tangential E and H at x = d and x - -d gives; 


Hg = HjjCo:h(vjjd-(|)g) 

(a) 

f!^ = HjjCosh(+(|)5) 

(b) 

V -- ^.^inh(-^s) 

S D 

(c) 

t 'i"h(V-*S> ' 

D m 

(d) 

(a) and (d) give 


r tanh(+v^d-*j) = 

D m 

(e) 

while (b) and (c) give 


V,. 

— — = tanh 4>c 

-Eg Vjj 

(f) 




(e) can be rewritten as 


ORIGINAL P;v2li 13 
OF POOR QUALITY 


2v^d - +4.^ + 4.5 


where we define by 


% = tann 


u-1 


m 


It is straightforward to evaluate the power propagated in the mode P 


P = }Re /; (fxH*).iidS 


= iE(,HbL^x(<t)xC 

where Ej^ and are the x and y component electric and magnetic field 
amplitudes respectively. C is the confinement factor given by 

m' m m S' s S' 


Re(w/3)/Ur 


where is the energy velocity = c/S, S being the slowing f|ct^r. 

For a symmetric structure 4 )^^^ *= and thus tanh(v^d) = -^ 77 \ or 
Eq. (III. 2. 8 ). ^ 


